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Abstract 


In  a  numerical  effort,  wall  jets  over  a  curved  wall  geometry  (Coanda  flows) 
have  been  investigated  using  Direct  Numerical  Simulation  (DNS)  and  a  new 
turbulence  modeling  approach.  The  present  research  is  motivated  by  exper¬ 
imental  results  indicating  that  the  presence  of  large  coherent  structures  can 
greatly  enhance  the  effectiveness  of  wall  jets  in  delaying  or  preventing ,  flow 
separation  on  airfoils.  For  utilizing  wall  jets  as  separation  control  devices  to 
their  full  potential  a  deeper  understanding  of  the  generation  and  development 
of  these  large  structures  is  essential.  The  objective  of  the  present  research  has 
been  to  investigate  the  effect  of  curvature  on  the  development  of  large  coherent 
structures  in  wall  jets,  in  particular  the  development  of  longitudinal  vortices 
(Gortler  type)  and  their  interaction  with  two-dimensional  vortical  structures. 
The  focus  has  been  on  Coanda  flows  over  circular  cylinders  using  two  differ¬ 
ent  computational  approaches.  In  the  first  approach,  the  new  Flow  Simulation 
Methodology  (FSM)  was  employed  on  a  body-fitted  grid  to  compute  a  turbu¬ 
lent  wall  jet  over  a  segment  of  the  Coanda  cylinder.  In  FSM,  the  contribution 
of  the  turbulence  model  is  ramped  in  depending  on  the  grid  resolution  rela¬ 
tive  to  the  local  Kolmogorov  length  scale.  The  turbulent  flow  was  generated 
through  a  bypass-type  transition  triggered  by  three-dimensional,  large  ampli¬ 
tude  forcing  of  a  laminar  wall  jet.  For  the  flat  plate  reference  case,  FSM  was 
employed  as  a  DNS,  LES  and  in  the  URANS  limit.  In  all  these  cases  FSM  was 
able  to  capture  the  large  two-dimensional  vortical  structures.  For  the  curved 
wall  geometry,  FSM  has  been  employed  thus  far  only  in  the  DNS  limit.  While 
Gortler  type  3D  vortices  are  observed  in  the  DNS,  they  remain  very  weak  com¬ 
pared  with  the  2D  structures.  This  was  probably  due  to  the  relatively  narrow 
computational  domain  chosen  for  the  simulations.  In  the  second  approach,  the 
complete  Coanda  flow  including  the  nozzle  was  computed  using  an  immersed 
boundary  technique  (IBT).  This  approach  includes  effects  of  the  nozzle  geom¬ 
etry  and  any  global  effects  of  the  flow  separation.  With  2D  simulations  the 
feasibility  of  IBT  for  computing  Coanda  flows  was  established.  Trends  for  the 
dependence  of  the  separation  location  on  the  Reynolds  number  measured  in 
experiments  were  confirmed  with  our  numerical  investigations. 
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1  Introduction 

Extensive  experimental  research  by  Wygnanski  and  co-workers  (Wygnanski  1997)  has 
provided  considerable  evidence  that  periodically  forced  wall  jets  (tangential  oscillatory 
blowing)  may  be  an  effective  method  for  the  control  of  separation  of  flows  over  single¬ 
element  or  segmented  airfoils  (including  flaps).  The  setup  of  two  such  experiments  is 
illustrated  in  Figure  1.  In  the  experiment  by  Seifert  et  al.  (1994),  a  high  momentum 
jet  that  is  blown  through  the  slot  along  the  upper  surface  of  the  flap  is  periodically 
forced  at  the  location  indicated  in  Figure  la.  The  additional  forcing  strongly  improves 
the  effectiveness  of  the  weak  wall  jet  in  preventing  separation,  and  thus  allows  for 
significantly  higher  angles  of  attack.  In  the  experiments  by  Seifert  et  al.  (1993)  on  an 
airfoil  with  a  straight  flap,  lift  increase  could  be  maximized  with  the  combination  of 
weak  steady  blowing  and  relatively  low  momentum  oscillations.  Coherent  structures 
generated  by  the  forcing  eliminated  the  large  wake  region  that  was  present  above 
the  wall  jet  when  only  steady  blowing  was  used  through  surface-normal  transport  of 
momentum. 

While  the  effectiveness  of  pulsed  blowing  for  separation  control  was  clearly  estab¬ 
lished  with  the  experiments,  many  of  the  fundamental  instability  mechanisms  that 
are  responsible  for  preventing  or  delaying  flow  separation  in  the  forced  flow  are  still 


a)  Experiment  by  Seifert  et  al.  (1994)  b)  NACA  0015  airfoil  with  slotted  flap 

from  experiments  by  Seifert  et  al.  (1993) 

Figure  1:  Experiments  on  active  separation  control  by  Wygnanski  and  coworkers. 
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not  understood.  For  example,  experiments  have  shown  that  wall  curvature,  both 
convex  and  concave,  as  well  as  tangential  blowing  (or  local  suction),  may  have  a  puz¬ 
zling  effect  on  the  effectiveness  of  periodic  excitation;  in  certain  instances  the  effect 
of  periodic  excitation  may  be  enhanced  while  in  other  situations  it  can  be  dimin¬ 
ished.  Understanding  the  behavior  of  the  large  coherent  structures  and  their  control 
for  flow  separation  on  airfoils  is  extremely  difficult  due  to  the  great  complexity  of  the 
complete  airfoil  flow  geometry  which  is  influenced  by  many  parameters;  e.g.,  adverse 
pressure  gradient,  curvature,  jet  momentum,  nozzle  geometry,  forcing  method.  For 
gaining  insight  into  the  major  contributing  instability  mechanisms  we  investigated 
prototypical  flow  geometries  that  include  some  but  not  all  aspects  of  the  airfoil  geom¬ 
etry.  In  previously  funded  research  (AFOSR  GRANT  F49620-97-1-0274)  the  effect  of 
curvature  was  excluded  by  focusing  on  pulsed  wall  jets  injected  into  separating  flat- 
plate  boundary  layers.  In  the  present  research  project  we  have  investigated  curvature 
effects  on  wall  jets  separately  by  studying  prototypical  wall  jets  over  curved  surfaces, 
specifically  Coanda  flows  where  fluid  is  blown  tangentially  along  a  cylinder  surface, 
and,  due  to  a  surface-normal  pressure  gradient,  adheres  to  the  curved  surface  over 
relatively  large  downstream  distances. 

As  illustrated  with  the  schematic  in  Figure  2  a  jet  is  issued  from  the  nozzle  at  angle 
&o  and  develops  along  the  cylinder  surface  as  a  so  called  Coanda  flow.  The  Coanda 


Figure  2:  Important  flow  parameters  for  the  Coanda  flow. 
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flow  adheres  to  the  cylinder  because  of  a  pressure  difference  between  the  ambient 
pressure,  Poo>  and  a  lower  pressure  at  the  surface,  ps,  which  builds  due  to  streamline 
curvature  and  can  be  computed  from  the  surface-normal  momentum  balance, 

r°°  u2  +  u'2  ,  .  . 

p~-p'=p/«  -RT7dr-  (1) 

In  the  experiments  of  a  turbulent  Coanda  flow  by  Neuendorf  and  Wygnanski  (1999) 
(experimental  setup  shown  in  Figure  3a),  the  jet  remains  attached  for  about  6  —  60  — 
210°.  Through  nondimensionalization,  the  surface  pressure  coefficient  is  defined  as 


cp  — 


(Poo  ~Ps)R 


...  —  f{6  —  #o,  Rex) , 


(2) 


(Po  ~Poo)b 

which  depends  on  the  distance  from  the  nozzle  and  a  so-called  “nominal  Reynolds 
number”  defined  by  Newman  (1961), 


Reft  — 


'(PO  -Poo)fr-R 

pu 2 


(3) 


With  the  following  definition  of  the  jet  exit  momentum  assuming  a  top  hat  profile  at 
the  nozzle, 


Jn  =  U2n  ■  b  «  —  (p0  -  poo) , 


26. 


(4) 


the  nominal  Reynolds  number  can  also  be  written  as 


R&N  ~  V  =  \^ReD  '  Reb  ’ 

which  is  the  geometric  mean  of  two  Reynolds  numbers  based  on  cylinder  diameter, 
D ,  and  on  nozzle  width,  6,  respectively. 

Due  to  streamline  curvature  the  Coanda  flow  can  be  susceptible  to  a  centrifugal 
instability  that  leads  to  Gortler-type  streamwise  vortices.  This  instability  depends 
on  another  modified  Reynolds  number,  the  Gortler  number 


Go  =  Re50.5^o.5/-R,  where  Reso.$  =  UmaxSo.5 jv 


(6) 


with  the  local  velocity  maximum  Umax  and  the  half  width  60.s  of  the  streamwise  ve¬ 
locity  profile  of  the  wall  jet.  In  the  experiments  by  Neuendorf  (2000)  (Figure  3a), 
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e  spacer  i  hot  wire  probe 


Figure  3:  From  experiments  by  Neuendorf  (2000)  on  the  Coanda  flow:  a)  experimen¬ 
tal  setup,  b)  contours  of  phase  averaged  streamwise  vorticity,  and  c)  profile  of  the 
spanwise  averaged  radial  velocity  at  9  =  80°  from  the  nozzle. 


it  was  convincingly  shown  that  longitudinal  vortices  of  the  Gortler  type  arise  in  the 
shear  layer  of  the  outer  flow  region  (Figure  3b, c)  by  virtue  of  a  centrifugal  instability 
mechanism  and  grow  in  size  and  strength  in  streamwise  direction.  This  3-D  mecha¬ 
nism  may  compete  and  interact  with  the  2-D  instability  mechanism  of  the  outer  shear 
layer  and  therefore  reduce  the  effect  of  any  2-D  periodic  forcing. 


The  main  focus  of  the  present  research  has  been  to  capture  and  investigate  the 
dynamics  of  the  2-D  and  3-D  large  coherent  structures.  Two  computational  avenues 
have  been  pursued  to  achieve  this  goal.  In  the  first  approach,  the  flow  is  computed  on 
a  body-fitted  grid  using  an  incompressible  Navier-Stokes  code  which  was  originally 
developed  for  Direct  Numerical  Simulation  (DNS)  of  transition  in  various  flows  over 
flat-plate  geometries  (boundary  layers,  plane  wall  jets)  and  which  has  been  extended 
to  curved  wall  geometries  (Zhang  and  Fasel  1999).  With  this  code  we  were  able 
to  perform  Direct  Numerical  Simulations  of  a  plane  turbulent  wall  jet  which  serves 
as  a  reference  (section  4.1)  and  of  the  flow  over  a  segment  of  the  Coanda  cylinder 
(section  4.2).  In  order  to  relax  the  resolution  requirements  for  computing  turbulent 
flows  a  new  turbulence  modeling  approach  has  been  developed  in  our  group  over  the 
past  few  years,  originally  in  cooperation  with  C.  Speziale  of  Boston  University  (see 
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Speziale  1998).  This  so  called  Flow  Simulation  Methodology  (FSM)  is  designed  to 
be  consistent  with  DNS  for  a  fine  grid  resolution  and  with  unsteady  RANS  for  a 
coarse  grid  resolution,  and  represents  an  LES  between  those  limits.  As  discussed  in 
section  4.3  we  could  demonstrate  the  capability  of  FSM  to  capture  large  coherent 
structures  in  wall  jet  flows  for  the  whole  range  of  application  from  DNS  over  LES 
to  URANS.  Simulation  of  the  whole  Coanda  flow  geometry  including  the  nozzle  and 
the  separated  flow  region  is,  however,  beyond  the  current  computational  approach 
on  a  body-fitted  grid.  Therefore,  a  Cartesian  grid  method  using  immersed  boundary 
techniques  (IBT)  for  representing  solid  bodies  inside  the  computational  domain  has 
been  employed  to  compute  the  whole  flow  field  around  a  Coanda  cylinder.  Results 
from  selected  test  calculations  of  this  novel  approach  are  presented  in  section  5.1  and 
a  discussion  of  results  for  the  Coanda  flow  is  provided  in  section  5.2. 
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2  Computational  Method  —  Body  Fitted  Grid 

The  Navier-Stokes  code  employed  for  our  numerical  investigations  on  the  body  fitted 
grid  was  originally  developed  in  our  research  group  for  boundary  layer  transition  (see 
Meitz  and  Fasel,  2000)  and  has  since  been  adapted  to  the  wall  jet  flow  geometry. 

2.1  Governing  Equations 

In  our  code,  the  complete  incompressible  Navier-Stokes  equations  are  solved  in  vorticity- 
velocity  formulation.  Regardless  of  which  form  of  time  averaging  or  spatial  filtering 
is  applied,  the  resulting  transport  equations  take  on  the  form, 

f  +  vx(sx9-v.f)  =  iv*3.  (7) 

The  overbar  on  the  velocity  vector  V  and  on  the  vorticity  vector  w  denotes  the  op¬ 
eration  of  averaging  the  appropriate  components  of  the  corresponding  vector.  For  sim¬ 
plicity,  overbars  are  omitted  from  now  on  in  this  report.  The  quantity  f  =  [r^]  is  the 
turbulent  stress  tensor  which  must  be  modeled.  Equation  (7)  is  non-dimensionalized 
by  the  free-stream  velocity  JJ0 0,  a  reference  length  L,  and  a  global  Reynolds  number 
defined  as  Re  =  Uo0L/u  where  v  is  the  kinematic  viscosity. 

For  the  present  calculations  of  the  Coanda  flow  a  simplified  curvilinear  grid  is 
used.  As  shown  in  Figure  4,  the  computational  domain  wraps  around  the  cylinder 
surface  with  a  constant  domain  height.  The  grid  coordinates  are  defined  as  follows: 

x  =  xo(s)  —  Tjsmd(s) ,  (8) 

V  =  y0(s)  +  r  cos  6{s) ,  (9) 

where  s  is  the  arc  length  on  the  surface,  xq (s)  and  yo(s)  are  the  surface  coordinates, 
and  r  is  the  normal  distance  from  the  surface.  This  simplified  curvilinear  grid  has 
properties  that  enables  us  to  retain  much  of  the  highly  efficient  and  accurate  nu¬ 
merical  scheme  employed  by  the  original  code  in  Cartesian  coordinates.  The  grid  is 
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Figure  4:  Schematic  of  the  computational  grid  for  Coanda  flow  computations  using 
a  body-fitted  grid  . 


orthogonal, 


dx  dx  dydy  _ 
dsdr  +  dsdr  ’ 


(10) 


and  the  resulting  metric  coefficients  are 
hi  =  h  = 


(i)  *(f) 


(11) 


where  k(s)  is  the  the  curvature  of  the  surface,  with  k(s)  >  0  for  concave  wall  and 
k(s)  <  0  for  convex  wall,  respectively.  Note,  /i2  =  /13  =  1. 


In  these  curvilinear  coordinates  the  vorticity  transport  equations  then  write  as 


dus  d  d 

~dt  +  dr^sV  ~UrU>  ~  q-(^zU-u}sw) 


duir 

~dt 


duz 

dt 


A 

i' 

A 

I®, 

hdsy 


2  2 k(s)  dur 


V'w,- 


h 2  ds 
_ld_ 
hds 


-  us 


72..  ,  2 K[s)dus 


Vwr  + 


OJt 


h?  ds  r 
v)  -  ~w<4 


r«2(s)i 

~Ur 

h2 

U)TU ) 

'k2(s)' 

h2  , 

w  —  UJZ 

■»)]  =  v 

r/c(s)/c'(s) 


h 3 


rK(s)n'(s) 


+ 


h2  Jj  ’ 


1  (12) 


h3 


+ 


f“]l  ■<-> 

(14) 
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where 

2  _  id_  (id\  i  d  (  o \  a2 

hds  \hds)  Jr  hdr  \  dr)  dz2 

rK'(s)  d  Id 2  /c(s)  5  d?_  _5^_ 

h 3  5s  +  h2  ds 2  h  dr  ^  dr 2  5z2 


The  vorticity  components  are  defined  as 

'du;  dv' 


U)< 


iOr 


/ 1  dv 

du 

OJz  =  - 

1  h  ds 

dr 

h  ds  ) 


+ 


uk(s)' 

h 


(15) 


(16) 

(17) 

(18) 


In  addition  to  the  vorticity  transport  equations,  the  following  Poisson  equations 


are  solved  to  obtain  the  velocity  field  in  curvilinear  coordinates: 


d2v  d2v  d2v 

ds2  dr2  dz2 

1  d  ( 1  du\  d2u 

hds  \hds  )  dz2 

1  d  ( 1  dw\  d2w 

hds  \h  ds  )  +  dz2 


d(hu>z)  du)s  d  ( 

ds  dz  +  dr  \ 

dujT  1  d  ( 1  dhv  \ 
dz  hds  \h  dr  ) 

1  duT  d  ( 1  dhv\ 
h  ds  dz  \h  dr  ) 


dhu 

ds 


1  du 
hds  + 


uk;(s)^ 

~r) 


(19) 

(20) 
(21) 


The  wall  vorticity  is  calculated  from  the  velocities  to  ensure  zero-divergence  of 


the  vorticity: 

d2cos  d2us 
ds 2  dz2 


LUT 

d(hujz) 

ds 


d2huT  d  ( d2v  d2v  d2v\ 
dsdr  +  dz  \5s2  +  dr2  dz2) 
d  d  f  dhu  l  du  vk(s)\ 

dz  dr  \  ds  h  ds  h  ) 

0, 

du/s  d  ( dhu  1  du  vk(s)\ 
dz  dr  \  ds  h  ds  h  ) 


/  d2v  d2v  32u\ 
\ds2  dr2  dz2  ) 


(23) 

.(24) 
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2.2  Flow  Simulation  Methodology  (FSM) 

In  our  turbulent  flow  calculations,  the  turbulent  stress  tensor  [r^]  in  equation  (7). 
is  modeled  using  a  new  Flow  Simulation  Methodology  (FSM).  FSM  has  been  de¬ 
veloped  in  collaboration  with  C.  Speziale  (1998)  as  a  turbulence  modeling  approach 
for  complex  turbulent  flows  containing  large  coherent  structures.  A  key  feature  of 
this  new  FSM  is  that  it  is  consistent  with  DNS  for  fine  grid  resolution  and  with 
Reynolds-Averaged  Navier-Stokes  (RANS)  calculations  for  coarse  grid  resolution. 
This  is  achieved  by  modeling  of  the  turbulent  stress  tensor  [t^]  in  equation  (7)  as 

Tij  =  f(A/Lk)r$ ,  (25) 

where  r-j  is  the  Reynolds  stress  tensor,  and  f(A/Lk)  is  the  so-called  contribution 
function.  The  terms  A  =  [( dx 2  +  dy2  +  dz2)/ 3]1/2  and  Lk  =  iW4/^/4  are  the  repre¬ 
sentative  computational  grid  size  and  the  Kolmogorov  length  scale,  respectively.  In 
the  coarse  grid  limit,  f(A/Lk)  1,  a  RANS  calculation  is  recovered.  In  the  fine 
grid  limit,  f(A/Lk)  — )•  0,  the  FSM  approaches  a  Direct  Numerical  Simulation  (DNS). 
For  the  Reynolds  stress,  r-j,  in  equation  (25),  the  Algebraic  Stress  Model  (ASM)  of 
Gatski  and  Speziale  (1993)  is  used  that  allows  for  a  misalignment  between  the  stress 
and  strain  rate  tensor  as  it  occurs  in  non-equilibrium  flows.  Further  details  of  the 
turbulence  modeling  in  our  FSM  approach  are  provided  in  a  paper1  presented  at  the 
AIAA  Aerospace  Sciences  Meeting  in  Reno  (Wernz  et  al.  2002). 

2.3  Numerical  Scheme 

A  detailed  review  of  the  numerical  schemes  used  for  the  solution  of  equation  (7)  and 
the  Poisson  equations  for  the  velocities  (in  Cartesian  coordinates)  can  be  found  in 
Meitz  and  Fasel  (2000).  With  the  simplified  curvilinear  coordinates,  the  numerical 
procedure  required  only  minor  extensions  by  including  a  few  extra  terms  into  the 


1 A  copy  of  the  paper  is  attached  at  the  end  of  this  report. 
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governing  equation.  In  summary,  the  time  dependent  spatially  varying  flow  field  is 
calculated  using  the  so-called  spatial  model  (inflow-outflow).  For  the  time  integration 
a  fourth-order  accurate  explicit  Runge-Kutta  scheme  is  employed.  Spatial  derivatives 
in  the  streamwise  direction  are  approximated  by  fourth-order  compact  differences 
employing  upwind  biased  and  downwind  biased  switching.  Wall-normal  derivatives 
are  computed  with  fourth-order  compact  differences  that  allow  for  grid  stretching. 
Pseudo-spectral  decomposition  is  utilized  in  the  spanwise  direction,  2.  Lastly,  a  fast 
Poisson  solver  (see  Meitz  and  Fasel  2000)  is  implemented  for  solving  the  Poisson 
equations  to  obtain  the  velocity  components. 

2.4  Boundary  Conditions  and  Disturbance  Generation 

For  calculations  on  the  body-fitted  grid,  the  nozzle  flow  is  not  computed.  Instead, 
a  laminar  wall  jet  is  prescribed  by  imposing  Dirichlet  conditions  at  the  inflow.  At 
the  wall,  no-slip  and  no-penetration  conditions  are  imposed.  The  upper  boundary 
is  assumed  to  lie  well  outside  the  wall  jet  where  the  flow  is  irrotational.  Hence  the 
vorticity  components  are  set  to  zero  along  this  boundary.  Also,  for  the  surface-normal 
velocity  a  Neumann  condition  is  prescribed  which  ensures  (through  continuity)  that 
the  surface-parallel  velocity  remains  zero  along  the  upper  boundary.  At  the  outflow 
boundary,  a  buffer  domain  is  employed  for  damping  out  the  time-dependent  fluctu¬ 
ations  of  the  flow  and  thus  minimizing  reflections  from  the  outflow  boundary.  The 
ramping'technique  is  similar  to  that  proposed  by  Kloker  et  al.  (1993)  and  is  discussed 
in  detail  by  Meitz  (2000).  In  spanwise  direction,  periodicity  is  assumed,  which  is  auto¬ 
matically  satisfied  by  the  spectral  treatment  of  this  direction  in  the  numerical  scheme. 

Close  to  the  inflow,  the  laminar  flow  is  tripped  with  a  volume  forcing  technique 
which  was  proven  very  effective  for  plane  wall  jets  (Wernz  2001).  Two  forcing  blobs 
are  placed  inside  the  inner  and  outer  region  of  the  wall  jet  at  two  wall-normal  locations 
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a)  distribution  of  volume  force  b)  time  signal  for  one  spanwise  Fourier  mode 


Figure  5:  Disturbance  generation  through  volume  forcing. 

where  the  flow  is  most  receptive  to  disturbances.  The  spatial  shape  of  these  forcing 
blobs  is  illustrated  in  Figure  5a,  3D  periodic  volume  forcing  (in  several  spanwise 
Fourier  modes)  at  large  amplitudes  (Figure  5b)  and  with  a  combination  of  different 
frequencies  for  each  spanwise  Fourier  mode  quickly  produces  a  broad  disturbance 
spectrum  and  leads  to  a  rapid  flow  transition  to  turbulence. 
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3  Computational  Method  —  Immersed  Boundary 
Technique 

The  second  computational  approach  is  a  Cartesian  grid  method  using  an  existing 
compressible  Navier-Stokes  code  developed  in  our  group  (Harris  1997).  This  approach 
provides  greater  flexibility  for  predicting  flows  with  complex  geometries  and  enabled 
us  for  the  present  research  to  model  the  Coanda  cylinder  from  the  experiments  by 
Neuendorf  and  Wygnanski  (1999)  including  the  nozzle  geometry.  To  account  for 
solid  bodies  inside  the  computational  domain  immersed  boundary  techniques  have 
been  incorporated  into  the  code.  These  techniques  employ  forcing  terms  which  are 
added  to  the  right-hand  side  of  the  momentum  equations  to  enforce  zero  velocity  at 
the  surface  of  a  body  or  throughout  a  solid  body. 


3.1  Governing  Equations 


The  governing  equations  for  the  compressible  code  by  Harris  (1997)  are  the  non- 
dimensionalized  compressible  Navier-Stokes  equations  in  conservation  form.  Nondi- 
mensionalization  of  the  governing  equations  is  carried  out  by  using  a  characteristic 
length  scale  (i.e.  the  cylinder  diameter,  D *2),  a  reference  velocity,  temperature,  den¬ 
sity,  and  specific  heat  (Z7*,  T0*,  p*0,  and  C*Q,  respectively).  Using  these  reference 
quantities,  non-dimensional  parameters  are  formed  that  characterize  the  compress¬ 
ible  flow,  namely  the  Mach  number,  Reynolds  number,  and  Prandtl  number, 
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Pr  =  (j.*qC;ok*q.  (26) 


The  non-dimensionalized  governing  equations  (continuity,  momentum  and  energy 
equations)  are: 
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2 The  asterisk  denotes  dimensional  variables. 
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+  [(j Et  +  p)Uj  -  UiTij  +  qj]  =  0  ,  (29) 

Ot  OXj 

where  the  total  Energy,  Et,  the  viscous  stress,  ri;-,  and  the  heat  flux,  qt,  are  defined 


(  T  Ukitk\ 

E‘  m  ',U-l)MS5  +  — ) 

p  ( du{  duj  2  duk 

Tlj  Re  \dx:  +  dxi  3  lj  dxk 


qi  (7-l  )Ma2RePrdxi' 

with  7  =  1.4,  and  Prandtl  number  Pr  =  0.71.  The  viscosity,  p,  is  calculated  using 


Sutherland’s  law, 


p(T)  =  T2/3 


1  +  5 
T  +  5 


where  5  = 


111  K 


Lastly,  the  equation  of  state  for  a  perfect  gas  is  invoked  to  determine  the  pressure, 


7  M  o2 


For  the  present  calculations,  the  two-dimensional  form  of  the  compressible  Navier- 
Stokes  equations  is  used.  In  vector  form,  they  write  as 

dU  dF  dG  TT  ,r 

at  ox  ay 

where 
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pv 

Q  =  PUV~Txy 

PV2+P-Tyy 

.  (Et  +  P)V  ~  UTxy  ~  VTyy  +  . 

The  right  hand-side  terms  in  equation  (35),  H  and  M,  represent  forcing  terms  which 
are  discussed  in  sections  3.2  and  3.3  respectively. 


3.2  Boundary  Conditions  —  Immersed  Boundary  Technique 

Boundary  conditions  for  the  intergration  domain  shown  in  Figure  6  include  a  small 
uniform  freestream  imposed  at  the  inflow  boundary  (left)  and  filtering  inside  a  buffer 
domain  at  the  outflow  boundary  (right)  to  damp  out  large  vortical  structures  before 
they  reach  the  boundary.  The  no-slip  and  no-penetration  conditions  along  the  sur¬ 
face  of  the  Coanda  cylinder  inside  the  computational  domain  are  satisfied  using  the 
immersed  boundary  technique.  Through  a  forcing  term,  H,  added  on  the  right-hand 
side  of  the  governing  equation  (35), 


Figure  6:  Schematic  of  the  computational  domain  for  Coanda  flow  computations 
using  a  Cartesian  grid  method. 
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the  fluid  particles  are  brought  to  rest  along  the  surface  of  the  cylinder.  For  computing 
the  components  of  H  several  methods  have  been  considered  including  a  feed  back 
control  method  proposed  by  Goldstein  et  al.  (1993)  and  a  direct  method  that  directly 
sets  the  desired  velocity  values  at  the  immersed  boundary  (Mohd-Yusof  1997)  have 
been  considered.  More  details  on  these  techniques  are  provided  in  the  paper  by  Terzi 
et  al.  (2001). 

3.3  Wall  Jet  Generation 

The  wall  jet  is  generated  by  placing  a  circular  mass  source  into  the  cavity  inside  the 
Coanda  cylinder  (Figure  6).  This  is  accomplished  via  the  second  forcing  term  on  the 
right-hand-side  of  governing  equation  (35), 

"  m  ' 

M  =  o  •  (37) 

_  h  J 

To  avoid  sharp  gradients  in  the  flow  field,  the  density  flux  distribution,  m,  is  repre¬ 
sented  by  a  Gaussian, 

rh(r)  =  A*— ^exp(ln/c(r/rm)2),  (38) 

*rm 

where  rm  is  the  radius  of  the  density  source,  r  the  distance  to  its  center  location, 
k  the  value  of  the  Gaussian  at  r  =  rm  (small),  and  Am  the  total  strength  of  the 
density  flux.  In  addition,  an  enthalpy  flux  h  =  mCvTm  must  be  included  in  term 
M  to  provide  the  added  mass  with  a  temperature,  Tm.  The  mass  source  causes  a 
pressure  buildup  inside  the  cylinder  cavity,  which  accellerates  the  fluid  through  the 
contraction  and  nozzle  exit  and  issues  the  wall  jet  along  the  outside  surface  of  the 
cylinder  (see  Figure  6). 
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Figure  7:  Example  of  a  computational  grid  for  Coanda  flow  computations  using  a 
Cartesian  grid  method. 

3.4  Numerical  Scheme 

The  compressible  Navier-Stokes  equations  are  solved  using  a  finite-difference  method. 
To  allow  for  comparison  with  incompressible  calculations,  a  Mach  number  of  Ma  = 
0.25  is  used.  The  equations  are  discretized  using  either  the  second-order  MacCor- 
mack  method  (Tannehill  1997)  or  a  fourth-order  split  finite-difference  discretization 
combined  with  a  fourth-order  Runge-Kutta  time  integration.  The  latter  method  is 
described  in  detail  by  Harris  (1997)  and  has  also  been  used  for  backward-facing  step 
flows  by  Terzi  (1998).  The  numerical  scheme  allows  for  clustering  of  grid  points  in 
both,  x  and  y  directions,  which  allows  for  grid  refinement  near  the  cylinder  inside  the 
computational  domain  (Figure  7). 
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4  Results  from  the  Body-Fitted  Grid  Calculations 

The  objective  of  our  Direct  Numerical  Simulations  of  the  turbulent  wall  jet  over  a 
curved  surface  is  to  capture  the  dynamics  of  the  2D  and  3D  coherent  structures  that 
develop  in  the  flow.  The  DNS  will  then  serve  as  a  reference  for  computationally  less 
expensive  FSM  calculations  that  compute  the  large  coherent  structures  while  model¬ 
ing  the  effect  of  small  scale  structures  with  the  turbulence  model.  In  the  calculations, 
a  laminar  wall  jet  with  Reso.5  —  16, 500  is  transitioned  to  turbulence  using  the  forcing 
method  discussed  in  section  2.4.  The  DNS  are  performed  on  an  equidistant  grid  with 
801  points  in  streamwise  direction,  on  a  stretched  grid  in  wall-normal  direction  with 
200  points,  and  with  21  symmetric  Fourier  modes  in  spanwise  direction.  For  the  FSM, 
the  resolution  is  reduced  in  the  streamwise  and  spanwise  directions. 

4.1  DNS  of  Turbulent  Plane  Wall  Jet  (Reference) 

As  a  first  step,  we  have  performed  a  DNS  of  a  turbulent  plane  wall  jet  with  the 
same  domain  size  and  resolution  as  used  later  on  for  the  Coanda  flow.  Due  to  the 
strong  spreading  of  the  turbulent  wall  jet  and  the  discrepancy  in  scale  between  the 
inner  and  the  outer  flow  region,  a  well  resolved  DNS  is  beyond  our  present  com¬ 
putational  capabilities.  The  current  coarse  DNS  therefore  represents  a  compromise 
using  a  fairly  narrow  domain,  a  symmetric  Fourier  transform  in  spanwise  direction 
(half  the  cost  of  the  full  transform),  and  a  coarse  wall-normal  resolution  in  the  outer 
flow  region.  Despite  these  shortcomings,  the  quality  of  the  computed  flow  is  surpris¬ 
ingly  good.  This  is  demonstrated  in  Figure  8  where  characteristic  time-averaged  flow 
quantities  from  the  calculation  are  plotted  in  inner  and  in  outer  coordinates  at  five 
streamwise  locations  and  are  compared  with  experimental  measurements  by  Eriksson 
et  al.  (1998).  Despite  the  short  time-period  for  the  averaging,  the  streamwise  and 
wall-normal  velocity  profiles  in  Figure  8a-d  are  quite  well  converged.  Except  for  the 
first  and  last  streamwise  station,  which  are  too  close  to  the  forcing  location  and  the 
buffer  domain  respectively,  both  streamwise  and  wall-normal  velocity  profiles  match 
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the  experimental  data  quite  well  in  both  inner  and  outer  coordinates.  The  agreement 
is,  in  fact,  far  better  than  in  any  previous  attempt  at  computing  a  fully  turbulent 
wall  jet  using  DNS,  particularly  for  y+  <  100  in  inner  coordinates.  Interestingly,  the 
streamwise  velocity  in  Figure  8a  exhibits  an  apparent  log-layer  which  is  much  more 
pronounced  and  closer  to  the  experimental  data  than  for  RANS  calculations.  The 
turbulent  shear  stress  in  Figures  8e&f  also  exhibits  a  very  good  match  between  DNS 
and  experiment,  particularly  close  to  the  wall  where  the  DNS  captures  the  near-wall 
asymptotic  behavior. 

The  presence  of  coherent  structures  in  the  turbulent  wall  jet  can  be  demonstrated 
with  the  snapshots  of  the  flow  field  presented  in  Figure  9.  Shown  in  Figure  9a  are 
color  contours  of  the  spanwise  average  of  the  spanwise  vorticity  for  the  computational 
domain  excluding  the  buffer  at  the  outflow,  shown  in  Figure  9b  for  the  same  domain 
are  iso-surfaces  of  two  levels  of  spanwise  vorticity,  u>z  =  —  1  (counterclockwise  rota¬ 
tion)  and  uz  =  1  (clockwise  rotation).  Due  to  the  3D  forcing  close  to  the  inflow, 
the  wall  jet  breaks  down  rapidly.  Despite  the  purely  3D  forcing,  large  2D  vortical 
structures  emerge  during  flow  transition.  The  structures  form  at  x  ~  20  as  dipolar 
vortex  pairs,  but  farther  downstream  only  the  vortices  with  counterclockwise  rotation 
in  the  outer  flow  region  persist,  while  the  vortices  with  clockwise  rotation  in  the  inner 
flow  region  have  mostly  disintegrated  into  small-scale  fluctuations.  This  clearly  indi¬ 
cates  that  the  outer  shear  layer  of  the  turbulent  flow  supports  2D  vortical  structures. 
The  development  in  time  of  these  structures  is  nicely  illustrated  with  the  space-time 
diagram  in  Figure  10a  which  contains  the  imprint  onto  the  wall  of  the  2D  structures 
inside  the  flow  field  over  time  (their  time-traces).  Shown  is  the  spanwise  averaged 
spanwise  wall-vorticity  normalized  by  its  local  time-average.  While  the  vortices  travel 
downstream,  subsequent  vortices  merge  and  form  larger  sized  vortices,  as  indicated 
by  the  merging  of  their  time- traces  in  the  space-time  diagram. 
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a)  streamwise  velocity  in  inner  coords. 


c)  wall-normal  velocity  in  inner  coords. 


y/8o.6 

b)  streamwise  velocity  in  outer  coords. 


e)  Reynolds  stress  in  inner  coords.  f)  Reynolds  stress  in  outer  coords. 

Figure  8:  Time-averaged  flow  from  a  DNS  of  a  turbulent  plane  wall  jet.  Shown  are 
the  streamwise  and  wall-normal  velocities  and  the  Reynolds  stress  in  inner  and  outer 
coordinates  at  five  streamwise  locations.  Also  shown  are  measurements  at  x  —  70 
from  experiments  by  Eriksson  et  al.  (1998). 
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Figure  9:  Snapshot  from  a  DNS  of  a  turbulent  plane  wall  jet.  Shown  is  the  instanta¬ 
neous  spanwise  vorticity  at  time  t=2.5/34Hz. 
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An  illustration  of  the  wall  vorticity  in  physical  space  is  provided  with  the  contour 
plots  in  Figure  10b  shown  at  three  time-instants  versus  streamwise  and  spanwise 
direction  (x-z  plots).  The  spanwise  direction  is  shown  to  scale  over  two  wave  length 
of  A^=1  =  2.  Since  the  flow  is  computed  with  a  symmetric  DNS,  the  wall  vorticity  is 
periodic  and  symmetric  in  spanwise  direction  causing  a  regularity  in  the  flow  pattern 
in  Figure  10b.  Close  to  the  forcing  location,  the  flow  is  highly  three-dimensional, 
but  2D  structures  begin  to  emerge  at  x  ~  12.5  and  repeatedly  double  in  size  as  they 
travel  downstream. 

The  streamwise  development  of  the  2D  vortical  disturbances  can  also  be  described 
with  a  Fourier  decomposition  in  time  of  the  spanwise  averaged  spanwise  wall  vortic¬ 
ity  for  each  streamwise  location.  In  the  two  contour  plots  in  Figures  11,  the  Fourier 
amplitude  of  the  wall  vorticity  and  of  the  normalized  wall  vorticity  are  plotted  as 
functions  of  frequency  and  streamwise  location.  In  Figures  11a,  the  spectrum  fills  up 
rapidly  downstream  of  a:  ~  12.  Over  the  interval  between  x  =  15  and  x  =  25,  high 
frequencies  dominate,  but  farther  downstream  a  clear  shift  in  the  spectrum  towards 
lower  frequencies  occurs  as  seen  in  Figure  lib. 

With  our  DNS  we  could  demonstrate  that  2D  coherent  structures  develop  in  the 
turbulent  plane  wall  jet  and,  as  the  jet  spreads  in  downstream  direction,  adjust  in 
size  through  vortex  merging. 
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a)  (x-t)  diagram  of  the  spanwise  average 
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b)  (x-z)  plots  at  three  time  instants 


Figure  10:  Evolution  of  coherent  flow  structures  in  the  turbulent  plane  wall  jet  gen¬ 
erated  through  bypass  transition.  Shown  is  the  spanwise  wall  vorticity. 
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Figure  11:  Fourier  spectrum  of  the  spanwise  wall  vorticity  versus  streamwise  direction. 
Shown  are  grey-scales  where  black  indicates  the  dominant  frequencies  in  the  spectrum. 
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4.2  DNS  of  Turbulent  Wall  Jet  over  Curved  Surface 

Encouraged  by  the  success  of  the  DNS  calculations  of  the  plane  turbulent  wall  jet  in 
section  4.1  we  have  performed  DNS  ofo  the  curved  wall  geometry  (Coanda  flow)  with 
an  otherwise  identical  computational  setup  with  the  same  domain  size,  resolution, 
inflow  conditions,  and  forcing  method  same  as  for  plane  wall  jet.  Only  small  steady 
volume  forcing  in  the  outer  flow  region  was  added  to  fix  the  location  of  any  longitu¬ 
dinal  vortices  that  might  develop.  The  radius  of  the  Coanda  cylinder,  R  =  0.8m,  in 
our  DNS  is  several  times  larger  than  in  the  experiments  by  Wygnanski  and  Neuen¬ 
dorf  (1999).  This  leads  to  a  Gortler  number3  of  Go  =  1300  at  the  inflow.  As  a 
result  of  choosing  a  large-sized  cylinder,  the  computational  domain  spans  only  about 
30°  of  the  cylinder  circumference  and  does  not  contain  the  region  of  flow  separation. 
The  present  computational  setup  was  chosen  due  to  numerical  limitations  of  our  cur¬ 
rent  curved- wall  Navier-Stokes  code  which  requires  a  constant  domain  height  and 
an  equidistant  grid  in  streamwise  direction.  Without  the  capability  of  adjusting  the 
domain  height  to  the  rapidly  spreading  wall  jet  the  computational  cost  of  calculating 
the  whole  Coanda  flow  is  too  great.  We  have  therefore  focussed  on  the  flow  region 
upstream  of  the  region  of  separation. 

The  Coanda  flow  computed  in  the  curved  wall  DNS  is  very  similar  to  the  plane  tur¬ 
bulent  wall  jet  discussed  in  section  4.1.  Snapshots  of  the  flow  field  at  one  time-instant 
are  provided  in  Figure  12,  with  color  contours  of  the  spanwise  averaged  spanwise  vor- 
ticity  shown  in  Figure  12a,  iso-surfaces  of  the  spanwise  vorticity  in  Figure  12b.  The 
2D  spanwise  vortices  in  the  outer  flow  region  of  the  wall  jet  are  clearly  discernible, 
their  dynamical  behavior  is  comparable  to  that  in  the  plane  wall  jet.  Again  the 
traveling  vortices  adjust  in  size  to  the  spreading  wall  jet  through  vortex  mergings. 

However,  due  to  curvature  the  Coanda  flow  is  also  Gortler-unstable  in  the  outer 
flow  region.  As  the  wall  jet  spreads  rapidly,  the  Gortler  number  does  increase  consid¬ 
erably  as  seen  from  the  plot  in  Figure  13.  Since  the  present  calculation  only  allows 


3  The  definition  is  given  in  equation  (6). 
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b)  iso-surfaces 

Figure  12:  Instantaneous  flow  field  from  the  DNS  of  a  turbulent  wall  jet  over  a  curved 
surface.  Shown  are  a)  color  contours  of  the  spanwise  average  of  the  spanwise  vorticity 
and  b)  iso-surfaces  of  spanwise  vorticity. 
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Figure  13:  Streamwise  development  of  the  Gortler  number  for  the  simulated  Coanda 
flow. 

for  symmetric  Fourier  modes  in  spanwise  direction,  a  meandering  of  Gortler-type 
streamwise  vortices  as  it  occurs  in  the  experiments  by  Neuendorf  (2000)  is  inhibited. 
We  therefore  were  able  to  isolate  the  streamwise  vortices  by  simply  time-averaging 
the  streamwise  vorticity.  Shown  in  Figure  14  are  the  color  contours  of  time-averaged 
streamwise  vorticity  at  six  streamwise  locations.  Also  plotted  for  reference  are  stream- 
wise  mean  flow  profiles  at  these  locations  and  the  half  width  of  the  wall  jet  which 
indicates  the  spreading  of  the  flow.  Distinct  vortices  exist  at  x  =  10  in  direct  re¬ 
sponse  to  the  forcing,  but  they  disappear  almost  completely  during  flow  transition 
(. x  =  20).  In  the  turbulent  flow  region,  however,  they  reorganize  and  are  clearly 
visible  at  x  =  30, 40.  Farther  downstream,  the  computational  domain  is  very  narrow 
compared  with  the  thickness  of  wall  jet  and  the  scale  of  the  vortical  structures  which 
therefore  cannot  increase  in  width  as  they  do  in  the  experiment.  To  perform  a  more 
realistic  calculation  the  domain  width  needs  to  be  doubled  or  quadrupled,  at  twice  or 
four  times  the  computational  cost.  Alternatively,  turbulence  modeling  can  reduce  res¬ 
olution  requirements  and  allow  for  efficient  calculations  with  a  wider  computational 
domain. 


r  [5mm] 
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Figure  14:  Development  of  streamwise  vortices  in  the  wall  jet  over  a  curved  surface. 
Shown  are  color  contours  of  time-averaged  streamwise  vorticity  at  six  streamwise 
locations,  the  streamwise  mean  velocity  profiles  at  these  locations,  and  the  half  width 
of  the  wall  jet. 
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4.3  FSM  of  Turbulent  Wall  Jet 

When  compared  with  DNS,  turbulence  modeling  using  the  Flow  Simulation  Method¬ 
ology  can  drastically  reduce  the  computational  effort  for  computing  turbulent  wall 
jets.  A  summary  of  the  results  from  our  FSM  calculations  will  now  be  given. 

The  DNS  from  sections  4.1  and  sec:BFcwj  represent  the  fine  grid  limit  of  FSM  for 
which  the  contribution  function  is  zero.  The  other,  the  coarse  grid  limit  is  an  Un¬ 
steady  RANS  where  the  contribution  function  is  one.  Many  features  of  the  turbulent 
wall  jet  can  be  captured  with  2D  URANS  calculations.  For  example,  in  a  turbulent 
wall  jet  a  mismatch  exists  between  the  location  of  the  maximum  of  the  streamwise 
mean  velocity  profile  and  the  zero  stress  location  which  is  closer  to  the  wall.  Con¬ 
sequently,  the  Boussinesq  approximation  does  not  hold,  and  for  modeling  this  effect 
using  conventional,  steady  RANS  complex  non-equilibrium  Reynolds  stress  models 
are  required.  As  discussed  in  the  paper  by  Wernz  et  al.  (2002),  the  URANS  calcula¬ 
tions  indicate  that  the  mismatch  between  stress  and  strain  results  from  a  distortion 
of  the  mean  flow  due  to  the  large  coherent  structures  (computed  in  URANS)  leading 
to  a  shift  away  from  the  wall  of  the  maximum  of  the  streamwise  mean  velocity  profile 
and  a  shift  towards  the  wall  of  the  zero  stress  location. 

URANS  also  provides  an  explanation  for  a  a  second  phenomenon,  the  repeated 
vortex  mergings  observed  in  DNS  of  turbulent  wall  jets.  Our  URANS  calculations 
of  the  forced  turbulent  wall  jet  suggest  that  subharmonic  resonance  occurs  between 
a  large  amplitude  primary  disturbance  and  a  small  amplitude  subharmonic  distur¬ 
bance.  In  the  URANS  calculation,  the  flow  is  forced  with  frequency  /  =  34 Hz  at  a 
large  amplitude  and  also  forcing  with  /  =  17Hz  and  /  =  8.5 Hz  at  a  much  smaller 
amplitude.  The  forcing  results  in  time-harmonic  propagating  disturbances  as  illus¬ 
trated  in  the  semi-logarithmic  plot  in  Figure  15.  Shown  for  the  wall  vorticity  are 
the  Fourier  amplitudes  of  the  fundamental  and  subharmonic  disturbances  versus  x. 
The  dashed  curves  represent  the  eigenbehavior  of  the  disturbances  which  exhibit  only 
moderate  growth  for  all  three  forcing  frequencies.  In  the  presence  of  the  large  ampli- 
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Figure  15:  Subharmonic  cascade  in  a  forced  turbulent  wall  jet  ( Rej  =  10, 000)  com¬ 
puted  with  URANS  (Wernz  2001).  Shown  versus  streamwise  direction  are  the  Fourier 
amplitudes  of  the  time  harmonics  represented  by  wall  vorticity. 

tude  primary  disturbance  with  /  =  34 Hz  (solid  blue  line),  however,  the  amplitude  of 
the  subharmonic  disturbance  (solid  red  line)  grows  rapidly,  and  surpasses  that  of  the 
primary  disturbance  leading  to  the  apparent  vortex  merging.  This  subharmonic,  now 
the  new  large  amplitude  primary  disturbance,  in  turn  amplifies  its  own  subharmonic 
(solid  green  line)  leading  to  a  second  merging  in  this  subharmonic  resonance  cascade. 

One  shortcoming  of  URANS  is  that  the  growth  rates  of  primary  disturbances  is 
lower  than  in  experiments.  This  is  likely  due  to  excessive  damping  of  time-dependent 
structures  by  the  turbulent  stress  from  the  turbulence  model.  In  3D  FSM  calculations 
where  the  value  of  the  contribution  function  adjusts  between  zero  and  one  depending 
on  the  local  grid  resolution,  the  turbulent  stress  term  contributes  less  than  in  the 
URANS  limit.  In  the  FSM  time-dependent  structures,  therefore,  grow  stronger  and 
reach  larger  amplitudes  than  in  the  URANS  calculations.  This  has  been  demonstrated 
with  a  case  study  for  FSM  presented  in  the  paper  by  Wernz  et  al.  (2002)  a  copy  of 
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which  is  attached  at  the  end  of  this  report.  In  this  case  study,  several  FSM  calculations 
with  various  grid  resolutions  and  consequently  different  model  contributions  have  been 
performed.  These  FSM  calculations  have  the  same  computational  setup  as  the  DNS 
in  section  4.1  with  the  same  3D  forcing  employed  to  trigger  flow  transition.  The 
turbulence  model  is  adjusted  once  breakdown  has  occurred.  Snapshots  of  the  flow 
field  from  five  FSM  cases  are  presented  in  Figure  16,  a  DNS  (case  FI),  an  LES- 
like  FSM  (case  F2),  an  FSM  with  a  coarse  (case  F3)  and  an  even  coarser  resolution 
(case  F4)  and  a  URANS  calculation  (case  F5).  Shown  are  color  contours  of  spanwise 
vorticity,  to  the  left  the  spanwise  averaged  vorticity  and  to  the  right  the  wall  vorticity. 
As  the  model  contribution  increases  from  case  FI  (0%)  to  case  F5  (100%),  the  small 
scale  structures  disappear  but,  except  for  the  URANS  in  case  F5,  the  large  scale 
structures  are  preserved. 

To  show  the  effect  of  the  turbulence  model  on  the  mean  flow,  in  Figure  17  the 
streamwise  mean  velocity  computed  with  the  various  FSM  cases  are  compared  with 
the  experimental  measurements  by  Eriksson  et  al.  (1998).  Shown  for  cases  F1-F5 
are  the  streamwise  velocity  profiles  in  inner  coordinates  at  one  streamwise  location. 
The  observed  trends  are  extremely  encouraging.  The  velocity  profiles  for  FSM  cases 
F2-F4  are  closer  to  the  experimental  data  that  the  profile  from  URANS,  case  F5.  In 
particular,  cases  F2-F4  exhibit  the  “log-like”  layer  which  is  measured  experimentally 
and  is  not  captured  by  conventional  RANS  calculations.  Case  F2,  the  LES-like  FSM 
which  has  the  same  grid  resolution  as  the  DNS  (case  FI),  appears  closest  to  the 


measurement  data. 
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Figure  16:  Instantaneous  flow  fields  from  the  FSM  calculations  (Cases  F1-F5)  at 
time  t  =  2/ZAHz.  Shown  are  color-contours  of  the  spanwise  average  of  the  spanwise 
vorticity  (left)  and  the  spanwise  vorticity  at  the  wall  (right). 
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Figure  17:  Mean  streamwise  velocity  profiles  from  FSM  calculations  (Cases  F1-F5) 
in  inner  coordinates.  Also  shown  are  measurements  from  experiments  by  Eriksson  et 
al.  (1998). 
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5  Results  from  Calculations  with  the  Immersed 
Boundary  Technique 

With  the  Cartesian  grid  method  using  immersed  boundary  techniques,  the  entire 
Coanda  flow  can  be  computed  in  one  calculation  including  the  nozzle  geometry  and 
the  region  of  flow  separation.  This  enables,  for  example,  investigations  of  the  effect 
of  the  nozzle  geometry  on  the  global  flow  field.  Traditionally,  immersed  boundary 
techniques  have  been  used  in  conjunction  with  computational  methods  of  low  order 
accuracy  (up  to  second  order).  Therefore,  a  strong  effort  is  being  made  in  our  research 
group  to  determine  the  impact  of  using  IBT  in  conjunction  with  our  higher  order 
methods  and  to  further  develop  IBT  so  that  overall  numerical  accuracy  is  preserved. 
Some  benchmark  calculations  are  now  presented  followed  by  a  discussion  of  our  results 
for  the  Coanda  flow. 

5.1  Validation  Cases 

Our  Cartesian  grid  approach  has  been  validated  with  three  key  benchmark  cases 
which  address  different  performance  features  of  the  method:  the  backward-facing 
step,  Tollmien-Schlichting  waves  in  a  flat-plate  boundary  layer,  and  the  cylinder  in 
a  free  stream  at  low  Reynolds  number.  Detailed  results  from  the  first  two  test  cases 
have  been  reported  in  the  paper  by  Terzi  et  al.  (2001). 

The  Backward-facing  step  is  a  standard  benchmark  case  for  CFD  codes  where 
the  reattachment  length  may  be  influenced  by  numerical  inaccuracies,  e.g.  excessive 
numerical  viscosity.  The  flow  physics  is  primarily  driven  by  the  shear  layer  that 
develops  from  the  step  corner  and  near-wall  accuracy  at  the  boundary  is  not  crucial 
for  obtaining  sufficiently  good  results.  We  have  compared  various  IBT  calculations 
using  different  boundary  treatments  with  the  flow  solution  from  a  previously  validated 
body-fitted  grid  calculation.  Results  from  two  of  these  calculations,  one  using  the 
feedback  control  method  by  Goldstein  et  al.  (1993)  the  other  the  direct  method 
by  Mohd-Yusof  (1997)  are  illustrated  in  Figure  18.  In  both  cases  (Figure  18b&c), 
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a)  Body-fitted  reference  case 
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x  X 

b)  Control  loop  (Goldstein  et  al.  1993)  c)  Direct  method  (Mohd-Yusof  1997) 


Figure  18:  Backward-facing  step  flow;  u-velocity  contours  for  body-fitted  grid  com¬ 
putation  and  for  two  IBT  methods,  immersed  boundary  techniques  applied  to  the  red 
area  only. 


the  IBT  is  applied  over  an  area  adjacent  to  the  step  base  marked  red.  Plotted  are 
contour  lines  of  streamwise  velocity,  u,  with  a  dashed  line  for  u  =  0  which  marks  the 
separation  point  at  the  wall.  Clearly,  both  IBT  results  are  in  excellent  agreement 
with  the  body-fitted  reference  case  in  Figure  18a.  In  particular,  their  reattachment 
points  match  up  very  well. 

In  contrast  to  the  backward-facing  step  flow,  near-wall  accuracy  is  crucial  for  cor¬ 
rectly  predicting  the  growth  of  Tollmien-Schlichting  (T-S)  waves.  Results  from  the 
calculations  by  Terzi  et  al.  (2001)  are  presented  in  Figure  19  with  an  illustration 
of  the  computational  setup  in  Figure  19a.  In  these  calculations,  various  IBT  meth¬ 
ods  with  and  without  wall  corrections  are  compared  with  linear  stability  theory  and 
a  DNS  on  a  body-fitted  grid.  While  the  wall-normal  amplitude  distribution  of  the 
streamwise  velocity  component  agrees  well  for  all  cases  (Figure  19b),  the  phase  differ¬ 
ence  between  the  streamwise  and  wall-normal  disturbance  velocity  components  close 
to  the  immersed  boundary  varies  strongly  for  different  methods  (Figure  19c). 
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a)  computational  setup 


LST-  linear  stability  theoiy 
B  -  body-fitted 
D  -  direct  method 
G  -  control-loop  method 
I  -  incompressible  code 
S  -  solid-body  correction 
R  -  reflection  correction 
C  -  one-sided  correction 


b)  u-disturbance  amplitude  distribution  c)  phase  difference  between  u  &  v  dist. 


Figure  19:  Tollmien-Schlichting  wave  simulation:  comparison  of  results  from  various 
IBT  methods  with  linear  stability  theory  and  results  from  DNS  with  body-fitted  grid. 
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Since  this  phase  difference  has  a  strong  influence  on  the  downstream  growth  of  the  TS- 
wave,  a  significant  spread  results  in  the  growth  prediction  from  the  different  methods 
(Figure  19d&e),  particularly  for  second  order  methods.  For  higher  grid  resolution 
(with  significantly  increased  computational  cost),  the  agreement  of  IBT  with  the 
DNS  using  a  body-fitted  grid  does  improve,  however.  The  conclusion  is  that  for 
computing  the  development  of  time-dependent  disturbances  that  are  driven  by  wall 
instabilities,  current  IBT  methods  are  less  accurate  than  higher-order  body-fitted  grid 
method  (used  for  the  calculations  in  section  4)  and  need  further  improvement. 

In  the  two  previous  tests  the  immersed  boundaries  are  aligned  with  the  grid,  which 
is  not  the  case  for  the  Coanda  cylinder.  We  have  therefore  tested  the  IBT  for  a  third 
standard  benchmark  case,  the  cylinder  in  a  uniform  freestream,  which  is  close  to  the 
Coanda  flow  geometry.  Using  either  IBT  variants  considered,  the  direct  method  and 
the  feedback  control  method,  a  Karman  vortex  street  is  obtained  (Figure  20).  The 
Strouhal  number  for  the  vortex  shedding  is  St  =  0.23  for  the  direct  method  and 
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spanwise  vorticity 


Figure  20:  Flow  around  cylinder  in  a  uniform  freestream  computed  using  IBT  (ReD  = 
400,  Ma  =  0.2).  Shown  are  color  contours  of  spanwise  vorticity  where  the  cylinder  is 
indicated  by  the  grey  shaded  area. 
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St  =  0.22  for  the  feedback  control  method  which  is  in  excellent  agreement  with  the 
literature  (e.g.  St  =  0.22  reported  by  Saiki  and  Biringen  1996).  This  leads  to  the 
conclusion  that  the  IBT  appears  capable  of  computing  separated  flows  over  curved 
surfaces  and  is  suited  for  computing  the  Coanda  flow. 

5.2  DNS  of  Coanda  Flows 

The  IBT  approach  has  enabled  us  to  compute  the  whole  Coanda  flow  field  including 
the  mass  supply  inside  the  cylinder.  In  our  calculations,  the  experimental  setup  from 
experiments  by  Wygnanski  and  coworkers  was  modeled  including  the  plenum  chamber 
inside  the  cylinder.  At  this  point,  2D  DNS  of  the  laminar  Coanda  flow  have  been 
performed.  In  the  present  simulations,  the  Coanda  flow  is  generated  by  applying 
a  mass  source  inside  the  plenum  chamber  of  the  cylinder.  Results  from  one  such 
IBT  calculation  are  plotted  in  Figure  21.  Shown  in  Figure  21a  are  color  contours  of 
spanwise  vorticity.  The  jet  exits  the  nozzle  at  6q  =  270°  and  remains  attached  up 
to  6sep  ~  150°.  As  shown  with  the  perspective  view  of  the  pressure  distribution  in 
Figure  21b,  the  mass  source  causes  a  pressure  build  up  inside  the  cylinder,  where  the 
pressure,  p0,  is  about  20%  above  ambient  pressure,  Poo-  Over  the  nozzle  contraction, 
the  pressure  drops  as  the  fluid  accellerates  and  the  surface  pressure,  ps,  is  below  p 
throughout  the  length  of  the  attached  flow  region.  The  Mach  number  (Figure  21c) 
peaks  at  Ma  =  0.4  at  the  nozzle  exit  but  then  drops  rapidly  below  0.3,  so  the  flow 
remains  mostly  incompressible.  Due  to  the  flow  acceleration,  the  temperature  inside 
the  cylinder  (Figure  21d)  drops  a  few  percent  below  ambient  temperature.  This  small 
variation  may  be  further  reduced  by  adding  the  mass  at  a  temperature  slightly  above 
ambient. 

In  a  parameter  study,  the  effect  of  the  nominal  Reynolds  number  on  the  length  of 
the  attached  flow  region  of  the  Coanda  flow  has  been  investigated.  The  goal  has  been 
to  keep  the  wall  jet  attached  for  as  long  as  possible.  Flow  parameters  for  the  cases 


Figure  21:  Laminar  Coanda  flow  computed  using  IBT.  Shown  is  the  stationary  flow 
field  from  case  3  (see  Table  1). 


are  listed  in  Table  l.4  For  case  1  (Figure  22a)  which  has  the  lowest  Reynolds  number, 
Res  =  75,  the  flow  spreads  rapidly  after  exiting  the  nozzle  and  separates  after  a  short 
distance  of  about  60°  from  the  nozzle.  In  this  case,  the  jet  is  not  energetic  enough  to 
remain  attached  for  very  long.  In  cases  2-5  (Figure  22b-e),  as  the  Reynolds  number 
is  further  and  further  increased,  the  separation  location  moves  farther  and  farther 
downstream.  From  our  experience,  and  not  surprisingly,  a  thin  jet  with  high  jet  exit 
momentum  works  best  in  keeping  attached.  In  fact,  in  inviscid  flow  calculations, 
where  Re^  —>  oo,  the  flow  adheres  to  the  cylinder  all  the  way  around,  for  360°. 

4 A  definition  for  pressure  coefficient,  Cp,  and  nominal  Reynolds  number,  Re\,  is  given  in  equa¬ 
tions  (2)  and  (3)  respectively. 
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case 

1 

2 

3 

4 

5 

6  (suction) 

7 

nozzle  width/radius,  b/R 

0.1 

0.1 

0.1 

0.05 

0.05 

in 

nozzle  angle,  do 

135° 

270° 

270° 

270° 

270° 

mil 

nominal  Reynolds  number,  Re ^ 

75 

250 

450 

740 

610 

610 

640 

pressure  coeff.,  Cp[d  —  do  =  80°) 

0.01 

0.07 

0.4 

0.13 

— 

separation  angle,  dsep  —  d0 

60° 

140° 

130° 

143° 

— 

Table  1:  Parameters  for  laminar  Coanda  flow  calculations  using  IBT. 


For  the  viscous  flow  computed  with  a  time-accurate  scheme,  however,  the  Coanda 
flow  becomes  more  and  more  unstable  to  two-dimensional  vortical  disturbances.  It 
becomes  increasingly  difficult  to  reach  a  stationary  state.  In  order  to  stabilize  the  flow, 
as  a  remedy,  wall  suction  has  been  applied  near  the  nozzle  exit  in  case  6  (Figure  22f) 
which  other  than  that  is  identical  to  case  5  (Figure  22e)  without  wall  suction.  As  a 
result,  the  amplitude  of  the  vortical  disturbances  is  drastically  reduced  and  the  flow 
remains  attached  farther  than  in  case  5.  How  destructive  the  two-dimensional  vortical 
structures  can  by  for  the  laminar  Coanda  flow  is  illustrated  for  case  7  ( Re jv  =  740) 
in  Figure  23,  where  two-dimensional  counter-rotating  vortex  pairs  repeatedly  disrupt 
the  flow  during  startup  and  propel  themselves  through  mutual  induction  away  from 
the  wall.  This  process  has  also  been  observed  for  laminar  plane  wall  jets  (Wernz 
and  Fasel  1996).  In  three-dimensional  calculations,  however,  this  2D  mechanism  is  in 
competition  with  3D  instabilities  in  the  flow  that  strongly  reduce  the  strength  of  the 
near-wall  vortex  and  thus  mostly  prevent  the  ejection  process  from  occurring  (Wernz 
and  Fasel  1999). 
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Figure  22:  Parameter  study  on  the  effect  of  the  nominal  Reynolds  number  Re^  oi 
the  Coanda  flow.  Shown  are  color  contours  of  spanwise  vorticity. 
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Figure  23:  Snapshot  of  the  flow  field  from  case  7  (see  Table  1)  where  no  stationary 
state  was  reached.  Shown  are  color  contours  of  spanwise  vorticity. 
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6  Summary  and  Outlook 

In  the  present  research  project  we  have  investigated  turbulent  wall  jets  on  Coanda 
cylinders  as  prototypes  for  separating  flows  over  curved-wall  geometries.  Two  com¬ 
putational  approaches  have  been  taken,  both  providing  valuable  insight  for  further 
research. 

The  flow  over  a  segment  of  the  Coanda  cylinder  has  been  investigated  using  Direct 
Numerical  Simulations  (DNS)  on  a  body-fitted  grid.  With  DNS,  the  development  of 
large  2D  coherent  structures  in  both,  plane  and  curved  wall  jets,  has  been  investigated. 
These  2D  vortical  structures  emerge  in  the  outer  flow  region  of  the  wall  jet  during 
flow  transition  and  persist  far  into  the  turbulent  flow  regime.  They  adjust  in  size  to 
the  spreading  wall  jet  through  repeated  vortex  mergings.  Using  URANS  calculations, 
evidence  was  found  that  subharmonic  resonance  may  be  the  mechanism  that  leads 
to  these  vortex  mergings.  For  the  curved  wall  jet,  longitudinal  vortical  structures 
have  been  found  in  the  outer  region  of  the  wall  jet.  These  Gortler-type  vortices 
also  grow  in  size  as  the  wall  jet  spreads.  Due  to  the  narrow  computational  in  the 
present  DNS,  however,  a  spread  in  spanwise  direction  is  inhibited,  which  possibly 
leads  to  unnaturally  narrow  vortices.  The  Flow  Simulation  Methodology  (FSM)  has 
been  adapted  for  wall  jet  calculations  as  a  new  approach  for  computing  this  non¬ 
equilibrium  flow  more  cost-effectively.  In  FSM  calculations  of  the  plane  wall  jet,  it  was 
demonstrated  that  the  correct  mean  flow  and  development  of  the  large  structures  can 
be  computed  with  coarsely  resolved  FSM  that  are  an  order  of  magnitude  less  costly 
than  DNS.  As  the  next  future  step,  the  FSM  now  will  be  applied  to  curved  wall  jets 
to  allow  for  wider  domains  in  spanwise  direction  than  in  the  present  DNS.  This  will 
allow  the  longitudinal  vortices  to  adjust  in  size  in  downstream  direction.  A  realistic 
simulation  of  the  flow  containing  both,  streamwise  and  spanwise  structures  will  enable 
us  to  study  their  interaction  and  actively  influence  their  role  in  the  flow  separation 
process  of  the  Coanda  wall  jet  farther  downstream.  Also,  two  current  computational 
limitations  of  our  Navier-Stokes  code  using  a  body-fitted  grid,  the  constant  domain 
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height  and  equidistant  grid  in  streamwise  direction,  must  be  overcome  to  employ  this 
approach  for  computing  the  actual  flow  separation  process. 

We  have  also  initiated  a  second,  novel  computational  approach,  a  Cartesian  grid 
method  using  the  immersed  boundary  technique  which  enables  us  to  compute  the 
entire  Coanda  flow  including  nozzle  and  separated  flow  region.  In  our  preliminary 
calculations  of  2D  laminar  flows  around  a  Coanda  cylinder,  we  were  able  to  keep  the 
flow  attached  up  to  about  140°  from  the  nozzle,  Major  trends  could  be  captured, 
namely  that  greater  jet  momentum  and  thinner  jets  lead  to  delayed  separation,  that 
larger  jet  momentum  favors  the  development  of  2D  coherent  structures,  which  in 
turn  promote  early  separation.  Since  near-wall  effects  may  have  an  impact  on  the 
flow  separation  process  further  efforts  will  be  made  to  improve  the  accuracy  of  the 
gradients  at  the  immersed  boundaries.  Also,  our  calculations  are  presently  limited  to 
2D  laminar  flow  due  to  computational  cost.  For  efficiently  applying  the  Cartesian  grid 
method  to  predict  flows  at  higher  Reynolds  numbers  and  turbulent  flows,  a  localized 
gridmesh  refinement  is  required  near  the  body  surface  and  other  flow  regions  with 
large  gradients,  e.g.  the  Adaptive  Mesh  Refinement  approach  (AMR).  Such  a  grid 
refinement  capability  is  currently  being  implemented  into  the  code. 
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Abstract 

The  evolution  of  large  scale,  coherent  structures  in 
turbulent  wall  jets  is  investigated  numerically  using 
our  Flow  Simulation  Methodology  (FSM).  Previous 
experimental  research  has  shown  that  these  struc¬ 
tures  significantly  change  the  mean  Reynolds  shear 
stress  distribution  and  also  the  skin  friction.  The 
combination  of  the  coherent  motion  in  the  shear 
layer  region  of  the  wall  jet  and  the  turbulent  bound¬ 
ary  layer  underneath  pose  a  unique  challenge  for 
turbulent  simulations.  We  use  FSM  to  investigate 
the  effect  of  these  large  coherent  structures,  both 
in  the  unsteady  Reynolds  Averaged  Navier- Stokes 
(URANS)  limit  and  in  the  Large  Eddy  Simulations 
(LES)  range.  In  addition,  Direct  Numerical  Simu¬ 
lations  (DNS)  will  be  used  to  show  the  emergence 
of  large  structures  with  strong  spanwise  coherence 
during  wall  jet  transition.  The  computational  re¬ 
sults  are  scrutinized  to  investigate  the  applicability 
and  performance  of  FSM  for  this  complex  turbulent 
flow. 

1  Introduction 

A  wall  jet  develops  when  a  fluid  is  blown  tangen¬ 
tially  along  a  surface.  Wall  jets  have  many  impor¬ 
tant  technical  applications  such  as  boundary  layer 
control  and  lift  enhancement  on  airfoils,  making 
use  of  their  tendency  to  adhere  to  curved  surfaces 
(Coanda  effect).  In  other  applications,  such  as  the 
cooling  of  turbine  blades  and  combustion  chamber 
walls,  a  wall  jet  is  used  to  cool  and  protect  surfaces 
from  hot  and/or  corrosive  fluids.  Despite  its  ap¬ 
plications  in  many  technical  systems,  the  wadi  jet, 
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especially  the  role  of  large  coherent  structures  in 
the  turbulent  wall  jet,  is  not  well  understood. 

With  respect  to  numerical  simulations,  the  tur¬ 
bulent  wall  jet  is  an  ideal  test  case  for  turbulence 
models,  since  it  combines  two  of  the  most  studied 
flows:  the  turbulent  wall  boundary  layer  and  the 
free  shear  layer.  However,  as  was  recognized  early 
on  [1],  the  interaction  between  these  two  prototypi¬ 
cal  flows  is  very  important  and  needs  to  be  modeled 
correctly  to  obtain  results  that  match  experimental 
findings. 

A  large  number  of  experimental  investigations  of 
the  mean  flow  of  turbulent  wall  jets  has  been  re¬ 
ported  in  the  literature  (see  review  articles  [2],  [1], 

[3] ).  Although  the  measurements  are  very  consis¬ 
tent  in  a  large  portion  of  the  mean  profiles,  the 
velocity  in  the  near  wall  region  is  extraordinarily 
difficult  to  measure  and  therefore  skin  friction  data 
vary  considerably.  The  most  important  parameters 
describing  a  wall  jet  are  its  local  velocity  maximum, 
Um ,  the  wall  normal  location  where  it  occurs,  ym, 
and  the  half  width,  t/o.5^  which  is  the  wall  normal 
distance  where  the  velocity  has  decayed  to  one  half 
the  local  maximum.  It  is  important  to  note  that 
in  the  turbulent  case,  ym  «  O.lSi/o.s,  and  the  edge 
of  the  jet  is  around  3r/o.s-  This  large  discrepancy  in 
length  scales  of  the  inner  boundary  layer  region  and 
the  outer  shear  layer  region  poses  unique  challenges, 
both  experimentally  and  numerically. 

Since  experimental  investigations  almost  exclu¬ 
sively  focused  on  the  mean  flow  quantities,  the  ques¬ 
tion  of  the  existence  and  influence  of  large  coherent 
structures  in  the  outer  shear  layer  on  the  mean  flow 
profiles  was  not  addressed.  From  a  hydrodynamic 
stability  point  of  view,  the  shear  layer  region  could 
be  unstable  to  disturbances  due  to  the  inflection 
point  in  the  mean  velocity  profile.  This  notion  was 
investigated  experimentally  by  Katz  and  coworkers 

[4] ,  who  have  shown  that  large  coherent  structures 
do  exist  in  the  turbulent  wall  jet,  and  although  these 
structures  are  visible  mainly  in  the  outer  region  of 
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the  wall  jet,  they  profoundly  change  the  whole  flow 
field,  as  indicated  by  the  significant  change  in  the 
mean  skin  friction.  More  importantly,  deliberately 
introducing  these  structures  using  periodic  forcing 
results  in  predictable  changes  of  the  mean  flow. 

While  there  is  an  abundance  of  experimental  data 
of  turbulent  wall  jets,  numerical  simulations  have 
only  recently  been  reported  ([5],  [6],  [7]).  The  dif¬ 
ficulties  in  computing  the  turbulent  wall  jet  can  be 
traced  back  to  the  discrepancy  of  length  scales  be¬ 
tween  the  boundary  layer  and  the  shear  layer  re¬ 
gions.  In  the  computations,  the  near  wall  resolu¬ 
tion  required  to  resolve  the  steep  gradients,  com¬ 
bined  with  the  large  domain  height  necessary  to 
accomodate  the  outer  shear  layer  part  makes  this 
flow  exceptionally  challenging  to  compute,  even  if  a 
turbulence  model  is  used. 

Most  numerical  investigations  to  date  have  fo¬ 
cused  on  the  mean  flow  quantities,  using  various 
turbulence  models.  The  results  presented  in  [4], 
however,  indicate  that  the  unsteady  coherent  mo¬ 
tion  strongly  influences  the  mean  flow  of  the  turbu¬ 
lent  wall  jet.  In  addition,  if  the  coherent  motion  is  a 
result  of  a  local  flow  instability,  then  steady  turbu¬ 
lent  simulations  are  unlikely  to  capture  the  relevant 
physics  correctly.  Therefore,  in  order  to  understand 
these  physical  mechanisms,  an  unsteady  simulation 
spanning  many  eddy  turnover  times  of  the  large  co¬ 
herent  motion  is  necessary;  at  the  same  time,  the 
simulation  has  to  be  accurate  enough  to  account  for 
the  random  turbulent  motion. 

To  overcome  the  computational  challenges  de¬ 
scribed  above,  various  techniques  have  been  devel¬ 
oped.  Extending  the  success  of  Direct  Numerical 
Simulation  (DNS)  from  the  laminar  regime  into  the 
turbulent  regime  proved  to  be  very  instructive  for 
understanding  the  relevant  physical  mechanisms  in 
prototypical  turbulent  flows  like  the  boundary  layer 
or  the  free  shear  layer.  However,  when  the  Reynolds 
number  is  increased,  these  simulations  quickly  be¬ 
come  too  expensive  even  for  todays  supercomput¬ 
ers. 

A  different  technique,  based  on  approximat¬ 
ing  the  turbulent  mean  flow,  is  Reynolds-averaged 
Navier-Stokes  (RANS)  simulation.  RANS  has  long 
been  used  in  industrial  applications  with  good  suc¬ 
cess.  Unfortunately,  RANS  cannot,  at  least  in  its 
original  form,  account  for  flow  unsteadiness  and  is 
therefore  limited  in  its  applicability  to  flows  where 


large,  coherent  (as  opposed  to  random)  motion  is 
not  important.  But  RANS  is  still  the  method  of 
choice  for  simulating  wall  bounded  shear  flows  for 
most  practical  applications. 

Large  Eddy  Simulation  (LES),  as  a  third  ap¬ 
proach,  is  able  to  compute  the  unsteady,  turbulent 
motion,  and  LES  has  become  the  method  of  choice 
in  research  for  flows  where  the  intricate  turbulent 
flows  physics  are  being  investigated.  Unfortunately, 
despite  extensive  research,  LES  could  not  live  up  to 
its  claim  to  be  a  method  capable  of  computing  all 
turbulent  flows.  For  a  discussion  of  the  shortcom¬ 
ings  of  LES,  see  ([8],  [9],  [10]). 

More  recently,  efforts  have  been  initiated  to  com¬ 
bine  the  strengths  of  these  methods  into  a  unified 
simulation  strategy,  sometimes  refered  to  as  hybrid 
models.  An  example  is  the  Detached  Eddy  Simu¬ 
lation  (DES)  [11].  In  principle,  DES  is  a  two  layer 
model,  combining  the  Spalart-Almaras  turbulence 
model  in  the  near  wall  region  with  a  traditional  LES 
away  from  the  wall.  The  switch-over  point  is  cho¬ 
sen  in  a  rather  ad  hoc  fashion,  which  is  not  critical 
if  large  coherent  motion  is  being  investigated  away 
from  the  wall  (like  in  wake  flow  simulations),  but 
unphysical  inconsistencies  can  be  introduced  if  the 
switch-over  point  is  not  chosen  carefully. 

Our  Flow  Simulation  Methodology  (FSM),  which 
we  have  been  developing  (originally  in  collaboration 
with  C.  Speziale  [9],  on  the  other  hand,  combines 
DNS,  LES,  and  RANS  in  a  consistent  manner  by 
using  the  same  underlying  turbulence  model  for  the 
RANS  and  the  LES  part  and  all  regimes  between 
RANS  and  LES  (e.g.  [12]).  In  fact,  FSM  main¬ 
tains  the  general  form  of  LES,  i.e.,  the  large  scales 
are  computed  and  the  small  scales  are  modeled. 
The  main  feature  of  our  FSM  is  that  the  subgrid- 
scale  stress  model  of  LES  is  replaced  by  the  prod¬ 
uct  of  a  contribution  function  and  a  state-of-the- 
art  algebraic  Reynolds  stress  model.  The  contribu¬ 
tion  function  is  designed  such  that,  depending  on 
a  ’physical’  numerical  resolution  (obtained  by  com¬ 
paring  the  grid  size  to  a  relevant  physical  length 
scale,  e.g.,  the  Kolmogorov  length  scale),  the  sim¬ 
ulations  approach  a  DNS  in  the  limit  of  fine  reso¬ 
lution  (or  low  Reynolds  number),  or  an  unsteady 
RANS  in  the  limit  of  coarse  resolution  (or  high 
Reynolds  number). 

In  this  paper,  our  Flow  Simulation  Methodol¬ 
ogy  (FSM)  is  used  to  investigate  the  origin  and 
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the  role  of  the  large  coherent  structures  in  wall 
jets.  FSM  has  been  developed  because  of  the  limi¬ 
tations  in  traditional  computational  strategies  such 
as  Direct  Numerical  Simulation  (DNS),  Reynolds- 
averaged  Navier-Stokes  (RANS)  computation,  and 
Large  Eddy  Simulation  (LES). 

2  Numerical  Method 


Note  that  the  contribution  function  depends  on  the 
ratio  (A/Lk),  i.e.,  the  physical  resolution  of  the 
computation  compared  to  a  length  scale  of  turbu¬ 
lence  in  the  flow.  It  should  be  pointed  out  that 
physically  relevant  length  scales  other  than  Lk  can 
be  chosen,  e.g.,  an  estimate  of  the  integral  length 
[14].  The  Reynolds  stress  t-j  in  equation  (2)  is  com¬ 
puted  using  the  anisotropic  Algebraic  Stress  Model 
(ASM)  of  Gatski  and  Speziale  [15], 


2.1  Governing  Equations 

The  governing  equations  are  the  three-dimensional, 
incompressible,  unsteady  Navier-Stokes  equations 
in  vorticity- velocity  formulation.  The  velocity  com¬ 
ponents  in  the  streamwise  (x),  wall  normal  (y),  and 
spanwise  (z)  directions  are  u ,  u,  and  u>,  respectively 
(figure  1). 

In  these  equations,  the  velocities  are  normalized 
by  the  jet  exit  velocity  Uj,  the  spatial  variables 
x,  y,  z  by  the  nozzle  width  5,  and  the  time  by  b/Uj. 
Taking  the  curl  of  the  momentum  equations,  which 
eliminates  the  pressure  gradient  terms,  yields  the 
transport  equations  for  the  vorticity, 

g  +  vx(ax?-v.?)-JLv’fl.  (i) 

For  simulations  using  a  turbulence  model,  f  = 
[rtj]  is  the  turbulent  stress  tensor.  In  the  two- 
dimensional  case,  this  set  of  equations  reduces  to 
one  equation  for  the  spanwise  vorticity  component 
ljz.  The  velocity  components  are  computed  from  a 
set  of  Poisson  equations  (see  [13]  for  more  details). 

In  order  to  close  equation  (1)  for  turbulent  simu¬ 
lations,  the  turbulent  stress  tensor  components 
need  to  be  modeled.  The  model  used  in  this  inves¬ 
tigation  is  our  new  FSM.  Following  Speziale  [9],  the 
turbulent  stress  tensor  components  are  written  as 

rl3.  =  /(A/L*)r«,  (2) 

where  is  the  Reynolds  stress  tensor  and  f(A/Lk) 
is  the  contribution  function, 

f(A/Lk)  =  [1  - 

exp(— 5  max(0,  A  —  2Lk)/NLk)]n-  (3) 

A  is  the  effective  computational  grid  size,  L ^  = 
(i^/e)1^4  is  the  Kolmogorov  length  scale,  and  N 
and  n  (in  this  study  n  =  1)  are  model  constants. 


\Kh 

-fin 7,0  ai-fSij 

+a2~  (SikWkj  +  SjkWki) 

~a3^2  (•S’ifcSfcj  - -W*)  .  (4) 


The  strain  rate  tensor  and  vorticity  tensor  are 
given  by 

■s.  _  I  fdui_  duA  W.._l(duL_dui 

tj  2  \<9xj  dxi )  7  tJ  2  \dxj  dxi 

o(5) 

The  function  /(t7,£)>  conjunction  with  k2/e , 
assumes  the  role  of  a  nonlinear  eddy  viscosity,  and 
is  given  by 

f(T]  0  = _ 1 _ w _ 

(6) 

where  the  expression  on  the  right  hand  side  of  equa¬ 
tion  6  is  a  regularized  form  of  f(r],  £),  which  is  used 
to  avoid  a  possible  division  by  zero  in  the  original 
expression  [9]. 

r]  and  f  are  invariants  of  the  ir rotational  strain 
rate  tensor  and  the  vorticity  tensor,  respectively, 


l£*3  k  rc  c  A1/2  c  a2k  /Tx7  TX7 

v  =  2^1(s'iS“]  =  ^77 {  “  v) 


The  values  of  model  constants  in  equations  (4) 
and  (7)  are  oq  =  0.227,  c*2  =  0.0423,  and  <23  = 
0.0396. 


The  turbulent  kinetic  energy  k  and  the  turbulent 
dissipation  rate  e  in  equation  (4)  are  computed  from 
the  standard  transport  equations  [16], 


Dk  dui  d  [  (vt  1  \  1  ,Q\ 

Dt  Ttj dxj  dxj  +  Re)  <9x; 
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Figure  1:  Computational  domain  with  schematic  of  volume  forcing  for  disturbance  generation. 


Dt  _  t  dui 
Dt  €l  kT%j  dxj 


Ce2fe2  J 


d 


The  eddy- viscosity  v?  and  the  wall  damping  func¬ 
tion  ft2  are  defined  as 


second-order  accurate  difference  stencils  are  em¬ 
ployed  for  the  streamwise  and  wall  normal  deriva¬ 
tives.  The  fourth-order  accurate  difference  method 
used  in  the  discretization  of  the  vorticity  transport 
equation  is  employed  for  the  spatial  derivatives  of 
the  turbulent  stress  tensor  components  (see  equa¬ 
tion  (1)). 


uT  =  CMfc2/€ ,  U  =  1  -  e-vflev'fc/io  (10) 

Finally,  the  constants  needed  for  the  solution  of 
equations  (8)  through  (10)  are, 

Cp  =  0.09 ,  Cel  =  1.44,  Ct  2  =  1.83, 

<rk  =  1,  =  1.3,  =  25. 

2.2  Numerical  Scheme 

A  detailed  description  of  the  numerical  schemes 
used  for  the  solution  of  equation  (1)  and  the  Pois¬ 
son  equations  for  the  velocities  is  given  in  [13].  For 
the  time  integration,  a  fourth-order  accurate  ex¬ 
plicit  Runge-Kutta  method  is  used.  The  deriva¬ 
tives  in  the  streamwise  direction  are  approximated 
by  fourth-order  compact  upwind  and  downwind  bi¬ 
ased  differences  at  consecutive  time  steps  and  wall- 
normal  derivatives  are  computed  with  fourth-order 
accurate  compact  differences.  In  the  spanwise  di¬ 
rection,  a  pseudospectral  approach  is  used.  Lastly, 
a  fast  Fourier  solver  is  implemented  for  the  solution 
of  the  velocity  Poisson  equations. 

In  addition  to  the  solution  of  the  vorticity  trans¬ 
port  equation  and  the  velocity  Poisson  equations, 
FSM  requires  the  simultaneous  integration  of  the  k 
and  e  equations.  For  these  equations,  a  second  or¬ 
der  ADI  method  is  used  for  the  time  integration  and 


2.3  Initial  and  Boundary  Conditions 

The  governing  equations  are  solved  in  the  compu¬ 
tational  domain  shown  in  figure  1.  At  the  inflow 
boundary,  Dirichlet  conditions  are  imposed  for  the 
vorticity,  the  velocity  components,  the  turbulent  ki¬ 
netic  energy,  and  the  turbulent  dissipation  rate.  For 
the  computations  presented  here,  the  inflow  condi¬ 
tions  are  obtained  from  a  solution  of  the  boundary 
layer  equations  for  the  turbulent  wall  jet.  At  the 
wall  y  =  0,  the  no-slip,  no- penetration  conditions 
are  imposed  on  the  velocity.  The  turbulent  kinetic 
energy  is  set  to  zero,  and  a  Neumann  condition  is 
used  for  the  turbulent  dissipation  rate.  At  the  free 
stream  boundary  y  —  ymaxi  &  decay  condition  is 
imposed  on  v,  and  the  wall  normal  derivative  is  set 
to  zero  for  all  other  quantities.  This  requires  the 
domain  height  to  be  large  enough  so  that  all  distur¬ 
bances  have  decayed  sufficiently  and  the  imposed 
boundary  conditions  do  not  influence  the  solution. 
Near  the  outflow  boundary,  a  damping  region  sim¬ 
ilar  to  the  one  proposed  by  [17]  is  used.  At  the 
outflow  boundary  itself,  all  second  derivatives  are 
set  to  zero. 

As  initial  condition,  the  vorticity  and  velocity 
components  are  prescribed.  The  turbulent  kinetic 
energy  k  and  the  turbulent  dissipation  rate  e  are 
prescribed  in  addition.  For  the  turbulent  simu¬ 
lations  (RANS  and  URANS),  these  quantities  are 
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obtained  from  the  solution  of  the  boundary  layer 
equations  for  the  turbulent  wall  jet,  otherwise,  the 
inflow  boundary  condition  and  the  initial  condition 
are  obtained  from  the  Glauert  similarity  solution. 

Disturbances  are  introduced  into  the  flow  using 
a  volume  forcing  technique  [21].  With  this  method, 
predominantly  vortical  disturbances  can  be  gener¬ 
ated  anywhere  in  the  flow  field.  In  the  cases  pre¬ 
sented  in  this  paper,  the  volume  forcing  is  applied 
near  the  inflow  boundary. 

3  Results 

3.1  Unsteady  Reynolds  Averaged 
Navier-Stokes  Simulations 

Two-dimensional  computations  of  the  mean  flow 
of  the  turbulent  wall  jet  are  used  to  validate  our 
new  FSM  in  the  limit  of  traditional  (steady)  RANS. 
These  computations  have  been  described  in  detail 
in  [7]  and  are  not  shown  here. 

Very  good  agreement  was  achieved  between  ex¬ 
perimental  findings  and  the  computational  results. 
However,  some  of  the  shortcomings  of  steady  turbu¬ 
lent  simulations  became  apparent.  First,  although 
anisotropic  effects  are  included  in  the  model,  the 
distribution  of  the  streamwise  Reynolds  stress  Rn 
showed  considerable  deviations  in  the  near  wall  re¬ 
gion.  Second,  we  found  that  the  non-equilibrium 
nature  of  the  turbulent  wall  jet  cannot  be  captured 
with  an  ASM  model  of  the  form  of  equation  (4)  (or 
any  other  equilibrium  k  —  e  model). 

3.1.1  Effect  of  forcing  on  the  turbulent  wall 
jet 

Using  the  Flow  Simulation  Methodology  (FSM)  in 
the  unsteady  RANS  limit,  the  effect  of  large  coher¬ 
ent  structures  on  the  turbulent  wall  jet  was  inves¬ 
tigated.  Even  though  FSM  was  developed  to  be 
used  as  a  subgrid  scale  turbulence  model,  it  is  of 
particular  importance  that  computations  in  the  un¬ 
steady  RANS  limit  yield  accurate  results.  The  dis¬ 
turbances  are  introduced  into  the  turbulent  base 
flow  using  the  volume  forcing  method  described  in 
section  2.3. 

For  comparison,  experiments  with  periodic  forc¬ 
ing  by  Katz  et  al  [4]  are  used.  They  achieve  a 
skin  friction  reduction  on  the  order  of  10%  with  a 


x/b=l00  x/b=1 50  x/b=2Q0 


Figure  2:  Effect  of  large  scale  structures  on 

the  mean  u-velocity  profile.  -  Base  Flow, 

•  ■  •  Case02, - Case03. 

forcing  amplitude  of  only  5%.  From  the  large  co¬ 
herent  structures  observed  in  a  flow  visualization, 
they  concluded  that  these  structures  result  from 
an  instability  of  the  turbulent  w'all  jet.  They  also 
showed  that  the  natural  unforced  turbulent  wall  jet 
has  a  strong  spanwise  coherence,  which  is  increased 
considerably  by  periodic  forcing.  If  these  coher¬ 
ent  structures  are  in  fact  nearly  two-dimensional, 
a  two-dimensional  simulation  should  be  sufficient 
to  model  the  effect  of  periodic  forcing  on  the  large 
coherent  structures,  and  as  a  consequency,  on  the 
mean  flow  of  the  turbulent  wall  jet. 

The  forcing  frequency  was  chosen  at  /  =  17  Hz. 
In  table  1,  the  forcing  amplitudes  for  the  cases 
described  in  the  following  are  summarized.  Fig¬ 
ure  2  compares  the  mean  flow  profiles  of  Case02 
and  Case03  with  the  profile  of  the  undisturbed  base 
flow  at  x/b  —  100, 150,200  (6  is  the  nozzle  height). 
In  Case03,  the  large  forcing  amplitudes  cause  a  sig¬ 
nificant  mean  flow  distortion.  Note  that  this  mean 
flow  distortion  increases  in  streamwise  direction,  in¬ 
dicating  that  the  disturbances  in  the  flow  actually 
enable  a  momentum  transfer  from  the  mean  flow  to 
the  large,  coherent  structures. 

To  illustrate  how  the  large,  coherent  structures 
influence  the  mean  flow,  the  Reynolds  shear  stress 
contribution  of  the  coherent  motion,  uVc,  the 
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Case 

CaseOl 

Case02 

Case03 

f/Uj 

5  x  1(T3 

1  x  1(T2 

2  x  IQ"2 

Table  1:  Disturbance  amplitudes  for  URANS  sim¬ 
ulations. 

Reynolds  stress  of  the  modeled  small  scales,  i?i2, 
and  the  sum  v!v' c  -f  R\2  for  Case03  are  plotted  in 
figure  3  as  a  function  of  the  wall  normal  distance. 
The  results  show  that  the  large  scale  structures  in¬ 
crease  the  Reynolds  stress  from  very  close  to  the 
wall  to  about  the  wall  jet  half  width  and  have  very 
little  effect  in  the  outer  part  of  the  shear  layer.  As 
can  be  seen  from  the  figure,  this  also  means  that 
the  y  location  of  the  zero  of  the  Reynolds  stress 
is  now  displaced  from  the  location  of  the  velocity 
maximum  towards  the  wall.  The  displacement  is 
in  surprisingly  good  agreement  with  experimental 
findings  (vr12= 0  ~  0.6y1/2  at  x/b  =  200).  This 
shift  has  so  far  been  attributed  to  nonequilibrium 
effects  and  its  origin  and  importance  has  been  the 
subject  of  considerable  controversy.  The  simula¬ 
tions  presented  here  indicate  that  the  large  un¬ 
steady  shear  layer  structures,  in  conjunction  with 
their  nonlinearly  created  harmonics,  are  a  major 
contributor  to  this  effect.  Figure  4  compares  R12 
from  the  ASM  model  without  large  scale  structures, 
the  sum  u'v'c  4-  R\2,  and  experimental  results  of 
Eriksson  et  al  [19]. 

In  summary,  ASM  in  the  limit  of  unsteady  RANS 
can  capture  the  effect  of  the  large  coherent  motion. 
However,  these  structures  have  to  be  introduced 
into  the  flow  by — in  this  case — periodic  blowing  and 
suction  because  unsteady  RANS  is  unsuitable  for  a 
computation  through  transition. 

3.2  Direct  Numerical  Simulations 

While  many  experimental  findings  are  matched  re¬ 
markably  well  by  our  URANS  calculations  (f=l  in 
eqn.  2),  some  considerable  differences  have  become 
apparent  as  well.  Firstly,  the  fully  developed  turbu¬ 
lent  flow  in  the  URANS  calculations  is  far  less  re¬ 
ceptive  to  forcing  than  the  transitional  region  close 
to  the  nozzle  where  forcing  occurs  in  experiments. 
As  a  result,  very  large  forcing  amplitudes  are  re¬ 
quired  for  producing  large  amplitude  disturbances 
and  numerical  instabilities  (due  to  the  turbulence 


x/b=  100  x/b=150  x/b=2C0 


Figure  3:  Reynolds  stresss  Ui2.  - random  part 

(modeled),  *  *  •  coherent  part, - sum  of  coherent 

and  random  parts. 


y 


Figure  4:  Reynolds  stresss  Ru  in  inner  coordinates. 

- Mean  flow  without  large  structures,  •  *  ■  Mean 

flow  with  large  structures,  o  experimental  results 
at  x/b  =  70,  •  at  x/b  =  150  [19]. 

model)  limit  the  achievable  disturbance  amplitude. 
Secondly,  in  URANS  calculations  the  generated  dis¬ 
turbances  are  mostly  damped  rather  than  amplified. 
This  is  in  contrast  to  the  experimental  results  by 
Katz  et  al.  [4],  who  reported  significant  disturbance 
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growth  in  the  fully  developed  turbulent  wall  jet.  On 
the  other  hand,  Schober  [20]  concluded  from  his  ex¬ 
periments  that  2D  coherent  structures  in  the  tur¬ 
bulent  wall  jet  are  damped  rather  than  amplified. 
In  his  view  these  structures  are  remnants  of  dis¬ 
turbances  originating  from  the  narrow  shear  layer 
upstream  at  the  upper  nozzle  lip,  where  the  flow  is 
still  transitional. 

To  address  these  issues  we  have  performed  Direct 
Numerical  Simulations  (DNS)  of  a  high  momentum 
wail  jet  through  transition  into  the  turbulent  regime 
(for  details  see  [21]).  In  these  simulations,  a  lami¬ 
nar  Glauert  wall  jet  (Re so. 5  =  16,500  at  the  inflow) 
is  perturbed  by  volume  forcing  with  several  (incom¬ 
mensurate)  frequencies  and  with  large  amplitudes. 
The  laminar  wall  jet  is  perturbed  close  to  the  inflow 
at  the  most  receptive  wall-normal  locations  in  both 
the  inner  and  the  outer  flow  regions.  Although  the 
forcing  is  applied  in  higher  spanwise  modes,  i.e.,  no 
two-dimensional  wave  component  is  introduced,  the 
simulations  show  the  emergence  of  large  structures 
with  a  very  strong  spanwise  coherence.  This  clearly 
indicates  that  the  wall  jet  naturally  supports  these 
structures.  From  a  hydrodynamic  stability  point  of 
view,  these  structures  could  be  connected  to  an  in¬ 
flectional  instability  of  the  outer  shear  layer  region 
of  the  wall  jet.  These  findings  are  in  good  qualita¬ 
tive  agreement  with  results  from  [20],  although  his 
experimental  setup  is  quite  different  from  our  sim¬ 
ulations.  A  snapshot  of  the  flow  field  is  shown  in 
figure  5,  where  the  isolevels  of  spanwise  vorticity  are 
shown  in  a  perspective  view  in  the  computational 
domain. 

To  determine  the  quality  of  the  turbulent  flow 
computed  with  the  DNS,  characteristic  time- 
averaged  flow  quantities  at  several  streamwise  loca¬ 
tions  have  been  compared  with  experimental  mea¬ 
surements  [19].  In  figure  6,  a  comparison  for  the 
streamwise  mean  velocity  and  the  Reynolds  shear 
stress  is  shown  in  inner  and  in  outer  coordinates. 
Except  for  the  first  and  last  streamwise  station, 
which  are  too  close  to  the  forcing  location  and  the 
buffer  domain  respectively,  the  velocity  and  shear 
stress  profiles  match  the  experimental  data  quite 
well.  The  agreement  is,  in  fact,  far  better  than  in 
any  previous  attempt  at  computing  a  fully  turbu¬ 
lent  wall  jet  using  DNS,  particularly  for  y+  <  100 
in  inner  coordinates.  Interestingly,  the  stream- 
wise  velocity  in  figure  6a  exhibits  an  apparent  log- 


layer  which  is  much  more  pronounced  and  closer 
to  the  experimental  data  than  for  RANS  calcula¬ 
tions.  However,  the  location  of  the  streamwise  ve¬ 
locity  maximum  in  the  DNS  is  farther  away  from 
the  wall  than  in  the  experiment.  The  reason  for 
this  might  be  that  the  2D  coherent  structures  are 
somewhat  more  intense  in  the  DNS  than  in  the  ex¬ 
periment  (and  in  our  URANS  calculations),  which, 
in  turn,  may  be  caused  by  insufficient  grid  resolu¬ 
tion,  particularly  in  the  outer  region  close  to  the 
edge  of  the  wall  jet. 

In  summary,  the  DNS  has  shown  that  coherent 
structures  are  generated  during  transitions  and  that 
these  structures  persist  in  the  turbulent  flow  regime. 
It  was  also  shown  that  the  mean  flow  data  from  the 
DNS  are  quite  close  to  turbulent  wall  jet  mean  flow 
data  obtained  from  experiments  and  RANS  calcu¬ 
lations.  The  agreement  is  especially  when  consid¬ 
ering  that  the  laminar  Glauert  wall  jet  used  as  in¬ 
flow  and  initial  condition  exhibits  a  vastly  different 
shape  and  streamwise  development.  The  goal  for 
our  FSM  investigation  is  to  explore  to  what  extend 
comparable  results  can  be  obtained  on  a  coarser 
computational  grid. 

3,3  Flow  Simulation  Methodology 

Extensive  investigations  of  the  FSM  have  been 
performed  for  a  zero  pressure  gradient,  flat  plate 
boundary  layer  [12]  and  good  agreement  with  the 
theoretical  predictions  and  other  (well  resolved) 
LES  and  DNS  simulations  was  found.  However,  for 
vrall  jets,  our  FSM  calculations  in  the  range  of  LES 
are  still  at  a  preliminary  stage.  We  have  chosen 
the  contribution  function  suggested  in  [12]  (eqn.  3). 
Here,  a  case  study  will  be  presented  to  demonstrate 
the  performance  of  FSM  throughout  its  range  from 
a  DNS  with  f(A,Lk)  =  0  (Case  FI)  to  URANS 
with  /(A,Ljt)  =  1  (Case  F5).  Case  FI  represents 
the  DNS  discussed  above  and  serves  as  a  reference. 
Cases  F2-F5  employ  the  DNS  as  an  initial  condi¬ 
tion  and  use  the  same  computational  setup  (com¬ 
putational  domain,  inflow  condition  and  forcing). 
However,  in  these  cases,  the  turbulence  model  is 
ramped  in  between  x  =  10  and  x  =  20  where  the 
flow  transitions  to  turbulence.  The  computational 
parameters  for  the  five  cases  are  summarized  in  ta¬ 
ble  2.  In  Case  F2,  which  has  the  identical  computa¬ 
tional  grid  as  the  DNS,  the  FSM  is  performed  in  the 


7 

American  Insitute  of  Aeronautics  and  Astronautics 


AIAA  2002-0589 


-volume  forcing 


filtered  COz-vorticity 
HHg§  C02=-1  (ccw) 


Figure  5:  Direct  Numerical  Simulation  of  wall  jet  transition.  Isolevels  of  spanwise  vorticity. 


regime  where  the  model  contribution  is  quite  small. 
Here,  the  turbulence  model  mainly  serves  to  reduce 
gridmesh  oscillations  that  are  prevalent  in  the  DNS 
as  the  resolution  in  the  DNS  is  not  quite  sufficient. 
Cases  F3  and  F4,  both  on  coarser  grids,  are  de¬ 
signed  to  have  virtually  all  small  scales  modeled  but 
without  the  strong  damping  effect  of  the  model  on 
the  larger  scales  as  in  the  URANS  limit.  Figure  7 
compares  the  level  of  the  model  contribution  for  the 
five  cases  at  one  time  instant,  t  =  2/34 Hz.  Shown 
are  the  values  of  the  contribution  functions  at  the 
wall  which  are  obtained  by  multiplying  by  an  ar¬ 
tificial  ramp  function  (e.g.,  see  figure  7,  Case  F5). 
Upstream  of  x  =  10,  the  model  is  shut  off  in  all 
cases. 

Results  from  the  case  study  are  presented  in  fig¬ 
ures  8-10.  An  impression  of  the  flow  fields  cor¬ 
responding  to  all  five  cases  is  provided  with  the 
isolevel  plots  in  figure  8,  which  display  the  span- 
wise  average  of  the  instantaneous  spanwise  vortic¬ 
ity.  For  Case  FI  (figure  8a),  the  flow  exhibits 
small  scale  structures  embedded  in  the  large  2D 
structures  throughout  the  computational  domain. 
The  large,  coherent  structures  consist  of  counter¬ 
rotating  vortex  pairs  corresponding  to  the  inner  and 


the  outer  flow  regions  of  the  wall  jet.  At  x  =  45, 
several  vortices  are  in  the  process  of  merging  and 
lifting  away  from  the  wall. 

For  Cases  F2-F4,  at  the  time  instant  shown  (t  = 
2/34 Hz),  the  flow  adjustment  following  the  intro¬ 
duction  of  the  turbulence  model  at  t  =  1  /MHz 
is  completed.  For  the  LES  in  figure  8b  the  smaller 
scales  are  reduced,  yet  the  large  scale  structures  still 
correlate  well  with  those  in  the  DNS.  The  structure 
at  x  —  45  is  somewhat  enhanced  and  in  the  process 
of  lifting  off  from  the  wall.  This  lift-off  of  dipolar 
vortices  by  mutual  induction  is  a  very  common  phe¬ 
nomenon  in  laminar  wall  jets  [22].  In  the  turbulent 
flow,  however,  the  near  wall  vortex  is  usually  too 
weak  to  propel  the  pair  away  from  the  wall  jet  core. 
With  the  contribution  of  the  turbulence  model,  the 
subgrid  scale  stress  replaces  the  weakening  effect  of 
the  small  scale  structures  on  the  large  structures. 
Indeed,  while  the  dipolar  structure  is  even  more 
pronounced  in  Case  F3  (figure  8c),  it  is  not  strong 
enough  to  lift  away  from  the  wall.  Interestingly,  the 
lift-off  process  does  not  occur  for  Case  F4  (figure 
8d),  which  has  a  very  coarse  spanwise  resolution. 
In  both  Cases  F3  and  F4,  the  smaller  structures 
are  mostly  confined  to  the  outer  edge  of  the  wall 
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a)  streamwise  velocity  in  inner  coordinates 


0  12  3 


y/50.s 

b)  streamwise  velocity  in  outer  coordinates 


c)  Reynolds  stress  in  inner  coordinates  d)  Reynolds  stress  in  outer  coordinates 


Figure  6:  Time-averaged  flow  from  a  DNS  of  bypass  transition  in  the  high  Reynolds  number  wall  jet. 
Shown  for  comparison  are  measurements  at  x  =  70  from  experiments  by  [19]. 


jet  where  the  contribution  function  shuts  down  the 
model  as  L k  increases.  Still,  size  and  spacing  of  the 
large  structures  is  approximately  preserved.  This  is 
clearly  not  the  case  for  the  URANS  in  Case  F5  (fig¬ 
ure  8e).  With  the  contribution  function  /  =  1,  the 
formation  of  the  large  2D  structures  from  3D  dis¬ 
turbances  is  suppressed  and  their  intensity  is  too 
small  to  be  visible  in  the  plot. 

The  increased  two-dimensionality  of  the  coherent 
structures  with  increased  model  contribution  is  il¬ 
lustrated  in  figure  9.  Shown  are  isolevels  of  spanwise 
vorticity  on  the  wall  (y  =  0).  For  Case  FI  (figure 
9a)  the  2D  structures  are  strongly  modulated  by 
3D  disturbances.  This  modulation  is  much  weaker 


for  Case  F2  (figure  9b)  where  the  small  scales  are 
reduced  in  amplitude  due  to  the  model.  For  cases 
F3  &  F4  (figures  9c&d),  the  modulation  is  hardly 
visible  and  for  Case  F5  (figure  9e)  all  structures 
are  suppressed  once  the  model  is  fully  turned  on 
(/  =  !)• 

The  temporal  development  of  the  2D  coherent 
structures  in  is  illustrated  in  the  space-time  dia¬ 
grams  of  spanwise  wall  vorticity  of  figure  10  for 
cases  F1-F4.  The  overall  pattern  of  the  time-traces 
for  the  vortical  structures  is  similar  in  all  cases,  i.e., 
the  spacing  (or  wave  length)  of  the  disturbances 
increases  in  streamwise  direction  as  the  wall  jet 
spreads.  The  most  pronounced  difference  is  the  wall 
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Case 

FI 

F2 

F3 

F4 

F5 

Type 

DNS 

LES 

coarse 

LES 

coarse 

LES 

URANS 

N 

oo 

2500 

1000 

1500 

0 

Ax[5mm] 

0.1 

0.1 

0.2 

0.2 

0.2 

Az(5mm] 

0.05 

0.05 

0.1 

0.25 

0.1 

i 

801 

801 

401 

401 

401 

k 

21 

21 

11 

5 

11 

Table  2:  Model  parameter  and  resolution  for  FSM 
cases. 

signature  of  the  lift-off  of  the  vortical  structures  in 
the  Cases  F2  &c  F3  (figures  9b  &c),  which  leaves 
wedge-like  voids  in  the  pattern. 

The  effect  for  Cases  F1-F5  of  the  turbulence 
model  on  the  mean  flow  is  illustrated  in  figure  11 
which  shows  the  streamwise  mean  velocity  in  inner 
coordinates.  The  observed  trends  are  extremely  en¬ 
couraging.  The  velocity  maximum  shifts  back  to¬ 
wards  the  experimental  measurement  curve  as  the 
model  contribution  is  increased  from  case  to  case. 
Cases  F2-F4  exhibit  the  “log-like”  layer  which  is 
measured  experimentally  and  cannot  be  captured 
by  conventional  HANS  calculations.  The  curves  for 
Cases  F3  FA  are  very  close  to  each  other  despite 
the  difference  in  spanwise  resolution.  Case  F2  ap¬ 
pears  closest  to  the  measurement  data.  It  is  unclear 
at  this  time  what  causes  Cases  F1-F4  containing 
time-dependent  structures  to  narrow  the  peak  at 
the  velocity  maximum.  Possibly,  the  narrow  width 
of  the  domain  and/or  the  fact  that  only  symmetric 
Fourier  modes  are  computed  might  play  a  role.  Al- 
ternaqtively  the  narrow  peak  may  be  an  artifact  of 
the  rapid,  bypass-type  transition  process  caused  by 
the  large  amplitude  forcing.  A  more  detailed  anal¬ 
ysis  of  the  present  results  and  further  investigations 
with  parameter  studies  will  be  needed. 

4  Concluding  Remarks 

In  this  paper,  the  development  of  large,  coherent 
structures  and  their  effect  on  the  mean  flow  of  a 
turbulent  wall  jet  has  been  investigated  with  our 
new  Flow  Simulation  Methodology  (FSM).  FSM 
combines  a  state-of-the-art  algebraic  stress  model 
(ASM)  with  a  contribution  function  /(A,  I*)  and 


Figure  7:  Wall  values  of  the  contribution  functions 
for  Cases  F1-F5. 

presents  a  unified  approach  to  computing  turbulent 
flows,  spanning  the  range  from  DNS  to  RANS.  The 
contribution  function  determines  the  appropriate 
amount  of  modeling  according  to  the  physical  reso¬ 
lution  of  the  simulation  by  comparing  the  local  grid 
spacing  with  a  local,  instantaneous  turbulent  length 
scale.  For  the  investigations  presented  here,  FSM 
was  used  in  the  limit  of  unsteady  RANS  (/  =  1), 
in  the  limit  of  DNS  (/  =  0),  and  for  two  levels  of 
modeling  contribution,  an  LES  and  a  coarse  LES. 

Without  adapting  the  underlying  ASM  model  or 
contribution  function  to  the  wall  jet  flow,  the  wall 
jet  mean  flow  as  well  as  the  evolution  of  the  large  co¬ 
herent  structures  observed  experimentally  are  cap¬ 
tured  using  FSM,  and  many  interesting  features  of 
this  complex  turbulent  flow  are  preserved.  Exam¬ 
ples  include  the  displacement  of  the  location  of  zero 
Reynolds  shear  stress  from  the  location  of  maxi¬ 
mum  velocity  (due  to  the  coherent  motion),  merg¬ 
ing  of  the  vortical  structures,  and  their  occasional 
lift-off  from  the  wall.  More  importantly,  these  intri¬ 
cate  physical  mechanisms  are  captured  with  a  sig¬ 
nificantly  reduced  resolution.  The  results  presented 
here  show  the  applicability  of  the  FSM  approach  for 
complex  turbulent  flow  simulations. 
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e)  Case  F5  (URANS) 


e)  Case  F5  (URANS) 


Figure  8:  Isolevels  of  instantaneous  spanwise  vor-  Figure  9:  Instantaneous  flow  fields  from  the  FSM 
ticity  (spanwise  average)  from  the  FSM  calculations  calculations  (Cases  F1-F5)  at  time  t  —  200/34 Hz. 
(Case  F1-F5)  at  time  instant  t  =  200 /MHz.  Shown  are  isolevels  of  spanwise  wall  vorticity. 
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a)  Case  FI  (DNS) 


spanw.average  <az  waN 


0  40  80 


b)  Case  F2  (“LES”) 


c)  Case  F3  (“coarse  LES”) 


d)  Case  F4  (“coarse  LES”) 


Figure  10:  Space-time  diagram  of  spanwise  wall  vorticity  from  the  FSM  calculations  (Cases  F1-F4)  at  time 
t  =  200/34 Hz. 
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Figure  11:  Mean  streamwise  velocity  profiles  from 
the  FSM  calculations  (Cases  F1-F5)  in  inner  coor¬ 
dinates  (x  =  40).  Shown  for  comparison  are  mea¬ 
surements  at  x  =  70  from  experiments  [19]. 
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